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ABSTRACT 

This learning module is concerned with the 
temperature fields the heat transfer rates^ and the coolant pressure 
drop in typical high temperature gas-cooled reactor (HTGH) fuel 
assemblies. As in all of the modules of this series , emphasis is 
placed on developing the theory and demonstrating its use with a 
simplified model. The heart of the module is the HTGH 
Thermal-Hydraulics Computer Code which solves for the radical 
temperature distributions in the fuel^ moderator^ and coolant at any 
axial station and then marches axially with an energy balance in the 
coolant. The code and its use are described in detail including a 
listing and definition of all variables, a discussion of all input 
requirements and resulting output, an annotated flow chart of the 
code, an explanation of all options in the code, and a listing of the 
code which gives enough comment statements to clearly indicate the 
operational steps being performed. By proper specification of the 
option, the code can either be used as an individual entity to study 
thermal-hydraulic aspects exclusively, or as a subroutine in the 
total HTGR module package to provide temperature feedback to the 
other modules. Examples are worked out using the code, (Author/SK) 
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THERMAL-HYDRAULICS ^K)DULE, TH-3 

HIGH TEMPERATURE GAS-COOLED REACTOR THERMAL-HYDHAULICS 

1.0 Object of Module 

The object of this module is to present the basic r. ' : . ■ of high 
temperature gas-cooled reactor (HTGR) thermal-hydraulics. requires 
the demonstration of ; 

• How the actual reactor core geometry can be modeled for 
simplified thermal -hydraulic analysis. 

• What information is necessary to characterize tV.^ thermal- 
hydraulic behavior of the reactor. 

• The development of the theoretical relations that permit 
the computation of these thermal -hydraulic characteristics. 

• The actual calculation of this information for the reactor. 
This calculation requires the use of the HTGR Thermal- 
Hydraulics Code, the description of which is included in 
this module . 

The thermal-hydraulic characteristics of the reactor are required 
for the determination of: 

• Fuel integrity 

• Moderator behavior 

• Coolant exit conditions 

• Helium compressor requirements 

• Temperature feedback for reactor neutronics calculations. 
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2,0 Content of the Module 

This learning module contains the thermal-hydraulics of high 

temperature gas-cooled reactors. Specifically, the module is concerned 

with the temperature field, the heat transfer rates and the coolant 

pressure drop in typical HTGR fuel assemblies. 

As in all of the modules of this series, emphasis is placed on 
developing the theory and demonstrating its use with a simplified model. 

The model is carefully selected to insure that analyses based on it will 
exhibit all of the important thermal-hydraulic trends of the typical reac- 
tor. The description of the core of a typical high temperature gas- 
cooled reactor and the modeling of its thermal-hydraulic characteristics 
are treated in the next section of this module. 

Following the geometry and modeling section, the basic theory 
governing the temperature distributions, heat transfer rates, pressure 
drops, and energy balance considerations is presented. The temperature 
distribution in the fuel is calculated assuming one-dimensional radial 
heat conduction. In the graphite moderator a two-dimensional finite 
difference calculation is used. The pressure drop in the coolant channels 
and the heat transfer coefficient for use in Newton's law of cooling are 
calculated from empirical relations developed for reactor coolant channel 
flows. Energy balances for small axial segments of the coolant channel 
are used to step the solution in the axial direction. Simple examples, 
illustrating the individual calculations, are worked out in detail for 
typical HTGR conditions. 

The heart of the module is the HTGR Thermal-Hydraulics Computer Code. 
Basically, the code solves for the radial temperature distributions in 
the fuel, moderator, and coolant at any axial station and then marches 
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axially with an energy balance in the coolant. The code and its use are 
described in detail. Included are a listing and definition of all vari- 
ables, a discussion of all input requirements and resulting output, an 
annotated flow chart of the code, an explanation of all options in the 
code, and a listing of the code which includes enough comment statements 
to clearly indicate the operational steps being performed. By proper 
specification of the options the code can either be used as an individual 
entity to study thermal-hydraulic aspects exclusively or as a subroutine 
in the total HTGR module package to provide temperature feedback to the 
other modules. Examples are worked out using the code. In typical 
examples, the location and magnitude of the maximum fuel temperature in 
the HTGR are found and the effect of undersized coolant flow channels on 
maximum fuel temperature and coolant outlet temperature are determined. 
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3.0 HTGR Fuel Geometry and Its Model 

The typical high temperature gas-cooled reactor core is approximately 
cylindrical In shape. A unit producing 3000 Mw(t) is composed of about 
4000 fuel element assemblies. These fuel element assemblies are in the 
shape of right hexagonal prisms about 31 in. high and about 14 in. across 
the flats. The fuel element assemblies are arranged in about 500 stacks, 
each 8 elements high, to form a core about 27 ft in diameter and 21 ft 
high. To provide a flattening of the neutron flux (and hence the power 
production) in the radial direction, the enrichment of the fuel in the 
fuel element assemblies is varied in cylindrical zones. Dowels are used 
to precisely align the fuel element assemblies. A small gap (about 0.04 
in.) exists between the elements at room temperature to accomodate thermal 
expansion. The core is surrounded by reflector and containment components 
which increase the inside reactor vessel size to about 38 ft in diameter 
and about 47 ft in height. 

The hexagonal fuel element assemblies are machined from graphite which 
acts as both the moderator for the neutrons and as structyral support for 
the reactor. The fuel, which generates heat internally during reactor 
operation, consists of highly enriched uranium and fertile thorium in 
carbide form. It is contained in fuel holes machined in the graphite. 
These fuel holes, which are typically 1/2 in. in diameter, extend vertically 
through about 95% of the 31 in. height of the right hexagonal prism. The 
cooling of the fuel element assemblies is provided by the downward flow of 
helium at about 50 atm. pressure through coolant holes machined between 
the fuel holes in the graphite. A typical coolant hole diameter is 5/8 in. 
The center-to-center spacing between fuel holes and between fuel and coolant 
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holes is typically 3/4 in. In each assembly there are about 200 fuel 

holes and about 100 coolant holes. The cross section of a HTGR fuel element 

assembly is shown in Figure 1. 

In each hexagonal fuel element assembly the six coolant channels in 
the corners of the hexagon are about 20% smaller in diameter than the other 
coolant channels in the assembly. Since these smaller channels amount to 
less than 6% of the total number, their smaller size will be neglected and 
the model coolant channel diameter will be selected as that of the central 
coolant channels, Tlie coolant flow conditions in the model channel will 
be taken as those for the average channel in the reactor. 

The uranium-thorium carbide fuel is coated with layers of pyrolytic 
carbon and bonded into rods. These rods fit into the fuel holes with about 
a 0,010 in, diametrical gap at room temperature. The rods span about 95% 
of the height of the fuel assembly element. To simplify the analysis, the 
model fuel rod is assumed to span the entire height of the element and the 
fuel rods are assumed to fit tightly into the fuel holes during reactor 
operation. 

The model fuel rod will have a variation in power density in the 
axial direction. This variation will be taken as 

(Z) = q'" cos (D 
o 

unless a computed actual axial variation is provided from another module. 
In Equation 1, q'" represents the thermal source strength per unit volume 
at any axial location Z, q'^' represents the thermal source strength per 
unit volume at the center of the fuel rod (Z=0) , and is the extrapolated 
height of the core. Both q'" and q'" are taken to be constant radially 
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Fifture 1. HTGR Fuel Geometry 
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throughout the fuel in the fuel rod. The magnitude of q'^' is- representa- 
tive of that of the average fuel rod in the core. The extrapolated height 
is calculated from 

H = H + 2 L (2) 
e e 

where H is the actual height of the reactor core and is an extrapolation 
length. The extrapolation length is the distance between the actual end 
of the core and the location where an extrapolation of the waveform re- 
presenting the actual neutron flux distribution within the core goes to 
zero. From neutron diffusion theory L can be shown to be about one 
migration length for a bare core. The migration length^'-can be calculated 
from neutron diffusion theory [1,2,3]. 

The triangular area shown hatched in Figure 1 represents the smallest 
area that can be analyzed to calculate the thermal-hydraulic behavior of 
the entire core. Note that all except about 1% (at the periphery of the 
hexagonal boundary) of the fuel element assembly can be constructed 
exactly from combinations of this basic "building block." The remainder 
can be approximated by this "building block." The segment of moderator 
enclosed in the triangular area will serve, in this module, as the model 
of the moderator for the thermal-hydraulic analysis. 

The reflector regions at the inlet (top) and outlet (bottom) of the 
core are each about 4 ft long. Therefore, for purposes of calculating 
coolant pressure losses, the additional lengths of coolant channel at the 
inlet and outlet of the core will each be assumed equal to 20% of the core 
height. 

*Numbers in brackets refer to items in References. 
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4.0 HTGR Thermal-HydgaiiljLc -Theory 

4.1 Internal Heat Generation 

As a result of nuclear fission in the fuel, heat is generated. The 
rate of energy generation in the fuel per unit volume is called the 
"volxjonetric thermal source strength," and can be calculated from 

q'*' - E^^^ 0^ ^ (3) 

where is the energy released per fission reaction (energy dimensions) , 

N is the fissionable fuel density (fissionable nuclei per unit volume) , 
f f 

is the effective fission microscopic cross section (dimensions of 
area), an(J ((> is the neutron flux (neutrons per unit area per unit time). 
Note that q'" has dimensions of energy per unit volume. The details of 
this calculation of volumetric thermal source strength are found in 
Reactor Statics Module 8. 

The volumetric thermal source strength varies throughout the reactor 
since ()>, and perhaps N^^, vary. For a cylindrical reactor, axial symmetry 
of the reactor fuel is generally a reasonable approximation. Therefore 
q'*' reduces to a function of only the radius, R, and the axial position, 
Z. Across any single fuel rod, the change in q'" is small (due to small 
change in R) and q"' can therefore be considered constant across its 
cross section. Thus, for any single fuel rod q*'* is a function of only 
the axial position. Also for any small axial segment of a fuel rod (the 
order of 2 to 3% of the total active length) the change in q''' is moderate. 
Using a constant average value for any such segment thus introduces little 
error into the analysis. This approximation is used in this module. 
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The heat generation in any small segment of the fuel can, be obtained 
by multiplying the volxjmetric theraal source strength by the volume of the 
segment. For steady state conditions this energy must be removed from 
the fuel. The mechanism by which this heat transfer occurs within the 
fuel is thermal conduction. 
Example 1 

A HTGR containing 5600 fuel elements stacked 8 elements high generates 
4000 l'fs.7(t). The cylindrical core of the reactor is 21 ft long and 32 ft 
in diameter. Each fuel element consists of a 3i.5 in. high right hexag- 
onal prism of graphite 14.2 in. across the flats and has 132 0.620 in. 
diameter fuel holes each containing a tight fitting uranium-thorium fuel 
rod. The volumetric thermal source strength of the fuel varies as q'" = 

q"* cos (ttZ/H ). The extrapolation length for the fuel is 1 ft. Find 
o e 

the average power per unit length (in kw/ft) and q'^' for the typical 21 ft 

high column of fuel rods in the core. 

Solution 

The power generated in the typical 21 ft high column of fuel rods is 

P total 4000 X 10^ /o onn /-J -J 1„ 

^ve = N =(5600) (1/8) (132) " 

. ave 43.3 kw o 
Thus, per unit length the average power is — — « 21 ft ' ^ ^'^^ ft* 

The power produced by the average fuel rod colxjmn can also be cal- 
culated by integrating the power produced in each differential volume of 
the fuel rods. Thus, 

H/2 

J q"' dV =J q'^' cos (ttZ/H^) (7TD^/4)dZ 



ave 



V -H/2 
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(H^ q';V4)l8in(irZ/H^)]^02 
(H q"*/2) 8in(irH/2H ). 



I » I 



Solving for q ^ 



q 



2P 

It. ave. 



° H sln(TrH/2H ) 
e r e 



(2X43,300 watt) (3.413 ^^^W^) 



(23 ft)(0.62. In)^ ^u^?^ sin[(ir)(21 ft)/(2)(23 ft)] 



4.86 X 10^ Btu/hr-ft^. 
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4.2 Temperature Distribution in the Fuel 

For conduction heat transfer the heat flux is proportional to the 
normal temperature gradient. V/hen the proportionality constant is inserted, 



where q is the heat transfer rate in the n-direction (energy / time) , 8T/8x^ 



area normal to the n-direction, and k is the thermal conductivity of the 
material (energy /degree -length-time) . Equation 4, which is Fourier's law of 
heat conduction, relates the heat transfer to the temperature field and is 
also the defining equation for the thermal conductivity. The thermal 
conductivity is a material property and its magnitude in general varies 
with the temperature of this material. Heat transfer properties of various 
reactor materials are tabulated as functions of temperature in References 
4, 5, 6, and 7. To reduce the complexity of heat transfer calculations, the 
thermal conductivity is often assumed to be constant and evaluated at an 
average temperature. The minus sign in Equation 4 assures that the heat 
transfer is in the direction of decreasing temperature. Equa^r.ion 4 shows 
that temperature gradients are required for heat transfer. In nuclear 
reactor applications, where there are high heat transfer rates, large 
temperature variations occur. One of the primary tasks of reactor thermal- 
hydraulic analysis is the prediction of this temperature distribution in 
the fuel . 

The heat conduction equation in cylindrical coordinates. 




(4) 



n 



n 



is the temperature gradient in the n-direction (degrees/ length) , is the 




q" ' ^ 1 ST 
k a 8T ' 



(5) 
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and the initial and boundary conditions prescribe the temperature distri- 
bution within the fuel rod. In this relation, the thermal conductivity, 
k, and the thermal diffusivity, a, have been assumed constant. The 
development of this equation is found in heat transfer texts [4, 5, 6] 
and several siiq)lified cases are given as exercises for the student. 
For steady state conditions, the unsteady term on the right hand side of 
the equation is zero, and for axial symmetry, the 6 variation disappears. 
For HTGR applications axial symmetry in the fuel rods is a reasonable 
approximation and since it greatly simplifies the analysis, this assumption 
will be made. It may also be observed that since the length of a fuel rod 
is much greater than its radius, the temperature gradients in the radial 
direction will be much greater than the temperature gradients in the axial 
direction. Therefore, to a good approximation, the heat transfer in the 
axial direction can be neglected with respect to that in the radial 
direction, and the resulting temperature distribution and heat transfer 
reduced to a one-dimensional case for any axial segment in which q'*' may 
be assumed constant. The governing differential equation for this case 
reduces to the ordinary differential equation. 

The solution of Equation 6 subject to the boundary conditions 

T » T at r " 0, 

(7) 

4^ « 0 at r = 0 

yields the temperature distribution in the fuel. The second boundary 
condition is obtained from the observation that the temperature distri- 
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bution must be continuous across the center of the cylinder. Observing 
that both q"' and k are constant and that 

Equation 6 can be written as 



dr2 



= .all! 

dr \ dr / 
Integrating twice results in 



r. 

f 



(9) 



T = - ^ # + C in r + r . ^^^^ 
4 1 2 

Applying the boundary conditions of Equation 7, the integration constants 

are 



= 0, 



2 o 



(11) 



Substituting into Equation 10 Rives 



T = T - r2. (12) 

o 4 k^ 

This relation shows the temperature distribution in the fuel to be 
parabolic with maximum temperature at the center. The heat transfer rate 
through any cylindrical shell can be calculated from Fourier's law which 
takes the form 

q - - k, A (13) 
^r f r dr 

where = 2TTr(AL) , AL being the length of the cylindrical shell. 
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Of particular interest are the temperature and heat transfer at the 
surface of the fuel; i.e., at r « r^ » 0^/2. At this location, 

T - T - -a-r- r/, (14) 
f o 4 f 



r - 



2 ^^'^ 
- irrj^(iiL)q"'. 

Noting that irr^ AL is the volume of the fuel rod segment of length AL, it 
is observed that the heat transfer out of the ^ » r^ cylindrical shell is 
indeed equal to the total energy generated as calculated from (q'*') 
(fuel volume). For one-dimensional heat transfer all of the energy 
generated within the fuel must be transferred out through the surface. 
Example 2 

Equation 14 relates the temperatures at the center and surface of a 
heat-generating cylindrical fuel rod. This equation contains the thermal 
conductivity of the fuel which, in general, is a function of temperature. 
However, the analysis leading to Equation 14 assumes the thermal conduc- 
tivity to be constant. Determine the temperature at which to evaluate 
to make the assumption of constant k^ compatible with a k^ that varies 
linearly with temperature. 
Solution 

The method of solution is to develop an expression relating T^ and 
T^ which assumes k^ = a + bT and then compare this result with Equation 
14. We begin with Fourier *8 law of heat conduction, 

1 A dT 
^ f r dr 



is 
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Recognizing that q = q'*' (volume) and = 2TTr(AL): 



q»'» TTr2(AL) = -k^ 2TTr(AL) ^ . 



Simplifying, separating variables, and substituting for k^: 

q' rdr = -2(a + bT)dT. 

, f f t 



Integrating from r = 0 to r = r^ with q***, a, b constant. 



It 



= [a + I (T^ + Xp] (T^ - 



q'" 



Comparing this with Equation 14 written as 



q'" vl 

T - T- = , . 
o f 4 



we see that the relations are equivalent for 



= a + |(T^ + Tp. 



te that this is precisely equal to k^ evaluated at the arithmetic 



No 

T + 

average fuel temperature; i.e., k(— 5_ )^ 
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Example 3 

In the UTGR Thermal-Hydraulics Code an expression for the average 
temperature of the fuel between the r = r^/2 and r = r^ radii is required* 
Show that this can be expressed in terms of the fuel centerline temperature, 
T^, and the fuel surface temperature, T^, as 

3 5 
f ave 8 o 8 f 

Solution 

The average fuel temperature in any region is defined as 



Expressing the differential volume and the temperature as functions of 
the radius, r, (using Equation 12 for the temperature) • 



T^ (T - 3^"^)27rrdr. 

f ave V I o ^ 

^r-/2 ^ 

2 

Noting that the volume contained between r^/2 and r^ is 3/4 Trr^ and 



integrating, 

■■I ^ l^f 

I r I 



i-nr^ L f J. 



f ave _^ 



77 b 



2 15^q' 

r. 



_ - ^ f 128 k 

5 q'"r/ 
T =— 

o 32 



But q" ' can be related to T^ and T^ by Equation 14 as 
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Therefore, 



T- = T 
r ave o 



8 o 
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4.3 Temperature Distribution In the Graphi te Moderator 

In theory, the temperature distribution in the graphite moderator 
can be obtained by solving the steady state heat conduction equation 
(without internal generation) subject to the boundary conditions prescribed 
along the periphery of the moderator. wever, the section of moderator 
of Interest is an Irregular shape as shown in Figure 1. Therefore, the 
bounding surfaces can not be expressed as lines of constant x, y, Q, r, or 
Z in either a Cartesian or cylindrical coordinate system and as a result, 
the boundary conditions are very difficult to apply. As a consequence, an 
analytical solution for the temperature distribution would be very difficult 
to obtain and even if obtained would be cumbersome to use. Thus, in the 
present module, finite-difference methods are used to determine the 
temperatures of Interest in the moderator. 

The basic finite-difference technique for determining the temperature 
distribution in a body consists of placing a nodal structure in the body, 
developing the system of algebraic equations that must be satisfied by the 
nodal temperatures, and solving the system of equations for the individual 

nodal temperatures. 

The nodal structure is selected as a compromise between extreme 
accuracy (many nodal points) and ease of solution (few nodal points) . 
The minimum number of nodes is restricted in the sense that there should 
be a nodal point at each location where the temperature is of interest. 
In general, nodal points are placed along the boundary of the region of 
Interest and then at equal spaclngs throughout the interior. However, as 
we will see later, the equal spacing is not a necessity, but simply a 
convenience. Each nodal point serves a finite-sized segment of the body 
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and has a well defined nodal volume and nodal surface area. The volume 
of all the segments taken together must comprise the body of interest. 
The temperature of the entire segment is taken to be equal to its nodal 
temperature. One convenient way of thinking of a node is as if all the 
material associated with each node is shrunk to the nodal point and the 
nodal point is connected to all adjacent nodal points by thin rods having 
the correct resistance to heat transfer. A numbering scheme must be used 
to provide identity for the individual nodes. Using the fact that there 
is no heat transfer through lines of symmetry, these lines can be replaced 
by adiabatic boundaries and the region that must be solved greatly re- 
duced. 

The algebraic equations satisfied by the nodal temperatures can be 
developed in two ways. The classic way is to recognize that the heat 
conduction equation must hold at the nodal point and then to express the 
partial derivatives that appear in the heat conduction equation as their 
finite-difference, approximation. Consider the nodal structure shown in 
Figure 2. It has been assumed that there are no changes in the Z" direction 
so that the problem reduces to a two-dimensional one. The shape of the 
body of interest suggested the use of the Cartesian coordinate system nodal 
structure that was set up. The setting of AX = AY was done for convenience. 
The heat conduction equation that must be satisfied at every point in the 
body, and therefore at node 9, is 

4^4^51!^ . 0. (16) 
3y^ k 

The partial derivatives in this equation are now approximated as follows: 
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Figure 2. Finite-difference Cartesian coordinate system example 
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AY 



^5 ^3 - ^ ^9 
Ay2 



(18) 



Substituting back Into Equation 16, 



0. 



AX 



AY' 



(19) 



For Ax = Ay, 



Tq + T,^ + + T,.. - 4 0. (20) 

8 10 5 IJ 9 k 



This is the desired nodal equation for node 9. Similar equations can 
readily be obtained for all interior nodes of the body. The extension 
to boundary nodes and especially convective boundary nodes such as nodes 
4 and 5 in this example is somewhat more difficult and will not be 
demonstrated for this method. Instead, the energy balance method, which 
is readily extended to all types of nodes, will be used in this module. 

To show the equivalence in the methods, the nodal equation for node 
9 will be determined by the energy balance method. For steady state 
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conditions the sum of the energy added to each nodal volume through its 
boundaries and the energy generated within the nodal volume must equal 
zero. Otherwise, there would be a net energy input to the nodal volume 
which violates the steady state (no changes with time) criterion. Thus, 
for node 9 , 

n + a + a „ + q„„ + q'" (nodal volume) = 0. (21) 

Each of the boundary heat transfer rates occur by conduction and can be 
determined from Fourier's law, 

q .-kA^-. 

n 9x 

n 

Expressed in finite-difference form this becomes 

^ ^ . 1, A ^ (22) 
n 

Applying this to each of the four boundary heat transfer rates, 

T - T 

-10^ = ^^^^ 

T - T 

= i'(^x)(i) 

qi3^ = k(AX)(l)^^V^ 
Note that the minus sign in Fourier's law has been absorbed in the AT 
and AX^ for each term. In general, one may note that the temperature of 
the node of interest always occurs behind the minus sign in the temperature 
difference when the heat transfer rates are into the node. This obser- 

26 
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vation saves much tedious chasing of minus signs in developing the nodal 
equations. Substituting the individual heat transfer rates into Equation 
21 and noting that the nodal volume is (AX) (AY) (1), 

^s'^^Q ^lo'^^Q ^5"^9 '^la-^g 

kAY + ^ +kAX +kAX +q' " (AX) (AY)=0 . (24) 

For AX = Ay this reduces to 

T +T +T+T -4T+ ^-T- (AX) = 0 (20) 
8 10 5 13 9 k 

which is identical to the result obtained earlier. 

The energy balance method will now be applied to node 6 to demonstrate 
how convective and insulated boundaries are treated. For node 6, 

^5->6 ^10-^ + q<^6 ^' " volrme) = 0. (25) 

For the two conduction heat transfer rates, 

T - T 

_ .AY. ""S ^6 



^5->6 ^ 2' AX 



T - T 
, ,AX, 10 6 



(26) 



'10-^ ^ 2' AY 
The convective heat transfer rate can be obtained from Newton's law of 
cooling, 

q = hA(T - T ) . (27) 

W 00 

For heat transfer into node 6, 

q^ = Mf)(T„- V. (28) 

Substituting these heat transfer rates back into Equation 25 and noting 



07 
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that the nodal volume is 

l^!^A+fcA|!l^A+h(f)(T.- V +q'''(^^ (29) 
For AX =■ AY, 

It +iT -(1 + ^) T, + . 0. (30) 

2 ^5 2 ^10 2k « ''^ ^ 2k ^ 6 4k 

This is the desired nodal equation for node 6. Note that the coefficient 

of T is the negative of the sum of the coefficients of the other tempera- 
6 

tures in the nodal equation. This observation can be used either as a 
check or to obtain the coefficient of the temperature of the node of 
interest. 

With the energy balanat ■•thod, the nodal equations for the 

remainder of the nodes in the example of Figure 2 can readily be obtained. 
With some practice, most of the nodal equations can be written by inspec- 
tion. Also, the extension of the method to three-dimensional and one- 
dimensional Cartesian coordinate problems should be obvious. 

Unfortunately, not all problems can be readily treated with Cartesian 
or even cylindrical coordinate system nodal structures. The present 
problem of determining the temperature distribution in the graphite 
moderator of the HTGR is such a problem. However, the energy balance 
method is quite general and can be extended to different shaped nodal 
volumes. Tetrahedron shaped nodal volumes are convenient for general 
three-dimensional problems and triangular right prism nodal volumes are 
convenient for two-dimensional problems since most bodies can readily be 
decomposed into a number of these shapes . Replacing arcs with their 
chords permits curved surfaces to be treated with little increased 

28 



25 

complexity and usually with acceptable accuracy. In the present moderator 
temperature distribution problem, temperature gradients in the vertical 
direction are several orders of magnitude less than those in the horizontal 
directions so that the problem can be treated as tv7o-diraensional. Thus, 
a triangular rif;ht prism nodal structure of unit depth (triangular in the 
plan view) will be used. 

An expression for the conduction heat transfer rates between nodes 
having a triangular structure must first be developed. Consider the 
general triangular nodal structure problem shown in Figure 3. The only 




T 

o 



Figure 3. Triangular nodal structure example. 
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restriction is that all of the angles of the triangles must be less than 
or equal to 90°. The nodal equation for node B is obtained from the 
following energy balance for steady state conditions: 



'Ia-^B 'IC^B 'Jp-^B Vb '^b'^ Sen ' °' 



(31) 



To calculate the internal energy generation for node B the nodal 

volume associated with node B must be defined. This nodal volume is 

taken as. the volume (per unit depth) enclosed by the dashed lines 

surrounding node B. This boundary is composed of perpendicular bisectors 

of all the lines linking node B to the other nodes . With the nodal 

volume of node B, V„, thus defined, the internally generated energy is 

6 

q -q'"V^. (32) 

The convective heat transfer rate can be calculated from Newton's 
law of cooling as 

q«^T. = h A (T« - T^) (33) 

^»^B conv B « D 

where A is the distance between the midpoint of the link between 

Conv B 

nodes A and B to the midpoint of the link between nodes B and C times 
the unit depth. 

To calculate the conduction heat transfer rates the concept of 
conduction shape factors is extremely useful. The conduction shape 
factor for heat transfer between nodes 1 and 2, S^^* defined as 

One sees immediately that for rectangular shaped nodal volumes the shape 
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factor is given by 

S„ = A /(AX ) (35) 
Rect ^ n' 

where is the heat transfer area perpendicular to the n-direction and 
AX^ is the spacing between the nodes in the n-direction. For triangular 
nodes it has been shown by Dusinberre [8] that the shape factors are 
related to the cotangents of the opposite angle of the triangles through 
which the heat transfer occurs. Specifically for path BG, 

he " 2 2 ^^^^ 

Similarly, assuming no conduction through the fluid on the convective 
boundary, 

he = I 

Writing similar expressions for the remaining conduction shape 
factors and substituting into Equation 31, 

^AB^V^B^ ^BC^V^B^ ^BF^^F-^B^ ^BG^^G-^B^ 

where 



= 1/2 cot BFA 

Ad 



S„- = 1/2 cot BGC 

S„_ = 1/2 cot BAF + 1/2 cot BGF 
BF 

S„^ = 1/2 cot BFG + 1/2 cot BCG. 
BO 



(39) 
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This relation can be simplified to 

T=.irST+ST+ST+ST + •'^""^ " T + ^-i ^1 (40) 

S„^^AB^A ^BC^C BF^F BG^G k » k ' ^ ' 



.Sionv B „ . J Bi 



B 

where 



hA „ 
conv B 

S„ = S_ + S„„ + S„„ + S„_ + 



'b " ''aB °BC ^BF " "BG k 
Similar expressions can be developed for each of the other nodes in the 
problem. 

The nodal equations developed above carry an implied assumption of 
constant thermal conductivity. If the temperature varies significantly 
from node to node and the thermal conductivity is a strong function of 
temperature the validity of this assumption is in question. In this case, 
a better thermal conductivity to use in each of the heat conduction 
terms is the arithmetic average of the thermal conductivities evaluated 
at the two nodal temperatures; i.e., 

^AB " ^^^V l^^V^^^- ^^^^ 

If k(T) is substituted into the nodal equations the resulting equations 
become nonlinear and more difficult to solve. However, often the use 
of a k evaluated at the average temperature of the problem suffices. 
Here a nodal volume weighted k would be the best value to use, but is 
tedious to calculate. Hence, a thermal conductivity evaluated at the 
arithmetic mean temperature in the problem is generally used. 

The nodal equations developed as described above for constant thermal 
conductivity constitute a set of linear algebraic equations for the nodal 
temperatures. These can be readily solved by standard, well-known methods 

32 



29 



such as matrix inversion, elimination, iteration, or relaxation. Each of 
these methods are well suited for computer calculations which are generally 
made for problems of more than a few nodes. 

Equation 40 is in a form directly useable for iteration methods. In 
iteration schemes, the temperature distribution in the body is first 
assumed. (Any rough, but logical, estimate works.) Then these assumed 
temperatures are substituted into the right-hand side (RHS) of the equations 
to calculate the temperatures on the left-hand side of the equations. 
Then the calculated temperatures are compared with the temperatures used 
in the RHS. If they agree to the desired accuracy, the problem is over; 
if not, the calculated temperatures are used in the RHS and the process 
is repeated until the desired accuracy is achieved. Two standard iter- 
ation schemes exist. In Jacobi (total step) iteration, the temperature 
in the RHS are updated after the entire set of equations has been used. 
In Gauss-Seidel iteration, the temperatures in the RHS are updated as 
soon as new values become available (from earlier equations in the set) . 
In general^ the convergence of Gauss-Seidel iteration is somewhat faster. 
Example 4 

A HTGR fuel element assembly has a cross section as shown in Figure 
1. The coolant hole diameter is D^, the fuel diameter is D^, and the 
spacing is S. At one level in the core the coolant temperature is T^, 
the heat transfer coefficient is h and the volumetric thermal source 
strength is q"'. Assuming no heat conduction in the axial direction, 
set up the equations for an iterative scheme to compute the temperature* 
distribution in the moderator. 
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Solution 

The section of moderator that must be analyzed Is the hatched area of 
Figure 1. An expanded view of this region with the nodal structure added 
Is shown In Figure 




Figure 4. Geometry of Example 4 showing nodal structure. 

The X and y coordinates of each of the nodal points are computed 
first In terms of the 3 geometric parameters, D^, D^, and S. 
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Y^ = 0 

Y2 = sinl 
Y3 = sin f 
Y^ = 0 

^■5 = ^2 

« 1 ^ TT 

Y^ = S cos - - 2 cosg- 
^8 = ^3 

Y5 = S cos f - I 



Tlie conduction shape factors are determined next. The cotangents 
in the shape factors are expressed in terms of the nodal point coordinates 
calculated above. 



S^2 4(Xi-X3)/Y3+lj^cotf 

m 

^14 " "2 [Arctan X^/Y^ + Arctan (X^ - 

^5=1 X3/Y3 4(X2-X^)/Y2 

m 

^25 " 2 [Arctan - X^)^^ + Arctan (X^ - X3)/Y3] 

+ J cot [Arctan (X2 - Xg)/(Yg - Y^) + Arctan (X^ - X3)/(Yg - Y3) ] 

^26 = i - ^5^/^^6 -^5^^i - V/<^3 " ^2^ 

^36 " 1 [Arctan^ (X^ - Xg)/(Y^ - Yg) + Arctan (X^ - X^)/(Y^ - Y^) ] 
+ J cot [Arctan (X2 - Xg)/(Yg - Y2) + Arctan (X^ - X2)/(Y2 - Y2) ] 
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S37 = i<^7- V/(^7- V 

= y cot [ir-Arctan Xg/(Yg - Y^) - Arctan /(Y^ - Yg) ] 
4(X3-X,)/(Y^-Y3) 
Sgg - I cot [Arctan Xj/CYg - Y^) + Arctan (X^ - " ^5) ] 

■^i^-^8^/^6 

S^g = I cot [f- - Arctan (Y^ - - Arctan (X3 - X^)/(Y^ - Y3) 

- Arctan (X^ - X^)/(Y^ - Y^) 1 
S^g = Y cot [u - Arctan (Yg - Y^) /X^ - Arctan (Y^ - " X^) ] 

^89 =iV(^9 - V- 

The nodal equations can now readily be developed using the above 
shape factors In the conduction terms. To assist In these developments, 
sketches of the more complex nodal volumes are useful In helping to 
visualize the heat transfer rates that must be taken into consideration. 
Use is also made of the observation that lines of symmetry can be treated 
as adiabatlc boundaries. 
Node 1: 

From the sketch for node 1 on the following page, 
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The last two terms may be combined by observing that they represent the 
energy generated in the entire 15° circular segment of the fuel. With 
this observation. 

Letting 

h " ^12 ^14 ^ h5 
and writing the nodal equation in a form convenient for iteration, 

1 '^^'"^f 

^1 - V^2^2 ^4^ ^ ^5^5 ^ -9rk- >• 
1 m 

Node 2: 

This node is similar to node 1 and its nodal equation can be written 
by inspection as 

^2 = bh2\ "-^ ^23^3 ^ ^25^5 ^ ^26^6 ^ ^TslT^ 
2 m 
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where 

^2 " ^12 ^23 hs ^26' 

Node 3: 

Again by inspection. 



where 



Node 4: 




Note that heat transfer between nodes 4 and 8 need not be considered 
because there is no finite area through which this heat transfer can 
occur. The use of image nodes (e.g. 5') across adiabatic boundaries 
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helps in visualizing the problem. Solving for T^, 



^ = i^'uh ^ ^45^5> 
4 



where 



h = ^14 ^5- 



Nodes 5, 6, and 8 are simple pure conduction nodes. 
Node 5: 



^5 = ?5^^15^1 ^25^2 ^45^ ^6^6 ^8^8^ 



where 



Node 6: 



^5 = ^15 ^25 ^5 ^58- 



^6 = ?/^26^2 V3 ^67^7 ^68^8 ^69 



where 



Node 8: 



where 



^8 ~ ^58 ^68 ^89\, 
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Node 7: 




Note that the actual circular arc heat transfer area is indicated. This 
will give a better representation of the heat transfer resistance in the 
convective layer than the chord length. 

V37(T3 - h'^^ Kh7^h - ^7> ^79(^9 - ^7>-^ ^ h ^c(^B " ^7^ = 0' 
Solving for T^, 

h ttD 

^7 = i(^37^3 ^ ^67^6 ^ ^79^9 ^ IaIT V 
7 m 

where 

h ttD 

^7 = ^7 ^ ^67 ^ ^79 ' 24ir ' 

in 

Node 9: 

h fTD 

^9 =VW6^ ^79^7 ■^^89^8-^ ''^^ 
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where 

h uD 

^9 = ^69 ^ '79 ^ ^89 ^ ITT • 

in 

This set of nine nodal equations can now be solved to determine the 
steady state temperatures at the nine nodal points in the moderator. 
This is the set of equations built into the HTGR Thermal-Hydraulics Code. 
In the HTGR Thermal-Hydraulics Code the equations are solved by an 
elimination method, rather than by iteration, to save computation time. 
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4.4 Temperature Drop in the Convective Layer 

All of the heat generated within the fuel ultimately must be 
transferred to the coolant. Thus, all of this energy must be transferred 
through the fluid layer near the surface of the moderator. This heat 
transfer mechanism, wherein the energy is carried away or convected from 
the solid surface by a fluid in motion is called convection. The heat 
transfer rate for convective heat transfer is related to the temperature 
difference between the surface and the bulk fluid, the driving potential 
for the heat transfer, by Newton's law of cooling, 

q = h A (T - T„) . (27) 
^c c c B 

This relation may also be taken as the defining equation for the heat 
transfer coefficient, h. The bulk temperature, T^ , is a mass-weighted 
average temperature of the fluid in the flow channel. It is formally 
defined by 



^pC^UT dA 



T = -J^-^- . (42) 

. B mC 

P 

This is the temperature that a thermometer would indicate if immersed in 
a cup of fluid collected from the discharge of the flow channel in the 
given location. 

Relations for calculating h from the flow characteristics and 
coolant properties have been developed. These will be treated in a later 
section of the module. For the moment, the quantity of interest is the 
temperature drop across the convective layer. In terms of h this can be 
obtained from Equation 27 as 



<?5 
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c B ttD (AL)h ^^-^ 
c 



where D is the coolant hole diameter, 
c 
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4.5 Pressure Drop in the Coolant 

All of the energy generated by nuclear fission in the fuel must be 
carried out of the reactor by the coolant which flows in the holes in 
the moderator. In a HTGR the coolant is gaseous helium which is force- 
circulated by a compressor. Enough pressure head must be provided by 
the compressor to overcome the pressure losses in the helium flow loop. 
The entire loop consists of the reactor flow channels, the heat exchangers, 
and connecting piping. Of interest here is the determination of the 
pressure loss incurred in the flow through the reactor core. TVo types 
of pressure losses will be considered. These are the frictional pressure 
loss along the coolant channel and the entrance and exit permanent 
pressure losses. 

The frictional pressure loss is a manifestation of the shear stress 
on the flowing fluid by the walls of the flow channel. This pressure loss 
is calculated from 



This relation may also be interpreted as the defining equation for the 
friction factor, f. The bulk velocity, U , found in Equation 44 is the 




(44) 



velocity averaged across the flow channel. It is defined by 




U dA. 



(45) 



This velocity is related to the mass flow rate by 



m = pAU, 



'B 



(46) 



41 



or to the mass velocity by 




A ^^B- 



(47) 



The value of the friction factor depends on the flow and surface 
conditions in the channel. The flow conditions are characterized by the 
Reynolds Number, 



It is well known that for Re below a critical value the flow is laminar, 
and for Re above the critical value the flow is turbulent. For internal 
flows this critical Re is about 2000. Laminar flow may be thought of as 
an ordered process in which fluid layers slide over one another, being 
retarded only by the molecular interaction between the layers. The vis- 
cosity of the fluid quantifies the magnitude of this interaction. In 
laminar flow any disturbance in the fluid is damped by the viscous 
action. In turbulent flow there is an additional random transport 
mechanism operable. This mechanism may be modeled as eddies (finite sized 
patches of fluid which retain their characteristics for finite times) 
moving throughout the fluid, transporting mass, momentum, and energy by 
virtue of their movement. This action is quite violent and results in 
transport rates much greater than those by purely - molecular activity in 
laminar flow. In turbulent flow, disturbances in the flow field grow 
and propagate resulting in increased turbulence levels downstream of the 
disturbances. Even in highly turbulent flows, there exists a layer near 
the solid bounding surfaces where the presence of the wall retards the 
penetration of eddies and therefore acts as a flow stabilizer. This 



(48) 



Re = 



^5 
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results in a laminar sublayer existing near any solid boundary. Even 
though this layer is very thin, much of the temperature and velocity 
change between wall and bulk conditions occurs in this layer. In reactor 
coolant channels the flow is highly turbulent with Re of the order of 
100,000 and greater quite common. 

The surface condition of the flow channel is characterized by the 
ratio of the surface roughness height to the diameter of the 

flow channel, e/D . The flow channel exhibits smooth tube behavior if 
c 

the roughness height is less than the thickness of the laminar sublayer4 
For reactor coolant flow channels the surface conditions are controlled 
so that this criterion is met. 

The smooth circular tube friction factor for the coolant channel 
can be determined from a ^toody chart where f is plotted versus Re . 
These charts are found in most fluid mechanics texts and handbooks [7, 
9, 10, 11]. An alternate method is to calculate f from a correlation in 
equation form. One of the most widely accepted for the turbulent flow 
conditions of reactor work is 

0.184 



circ „ 0.2 
Re 



Tliis relation, like the Moody diagram, was obtained from a curve fit to 

experimental data in smooth circular tubes. The friction factor defined 

in this module is the Darcy-Weisbach friction factor. Care must be taken 

not to confuse this friction factor with the Fanning friction factor 

(f = f , /4) also found in the literature. 
Fan D-W 

Tlie permanent pressure losses at the inlet and exit are the result 
of increased viscous energy dissipation resulting from increased turbulent 

i6 
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activity. These losses may be calculated from 

c 

In this relation K is the resistance coefficient for the inlet or exit 
geometry. The inlet may be approximated as a sudden contraction from a 
very large diameter to the coolant hole diameter. For this sudden 
contraction, the resistance coefficient is 0.5 [11]. Similarly, 
the exit may be approximated as a sudden expansion from to a very 
large diameter. For this expansion, K is 1.0 [11]. 
Example 5 * 

The HTGR described in Example 1 has an inlet coolant temperature of 
600 F, an inlet coolant pressure of 701 psia, and a coolant mass velocity 
of 8.3 X 10^ Ibm/hr-ft^ through coolant holes 0.826 in. in diameter. The 
coolant exit temperature is 1400 F. Estimate the frictional pressure loss 
in the core. Neglect entrance and exit losses. 
Solution 

To obtain a one-step estimate, the coolant properties are evaluated 
at the average coolant temperature and pressure. The average coolant 
temperature is 1000 F. The average pressure will be assumed to be 700 
psia. (This must be checked later.) At these conditions the helium 
properties required in the solution are 

p « 0.0889 Ibm/hr-ft 
P (700 ^) (144 

(386 ^^)( 1460 R) ft"* 
The average coolant velocity is found from 

\j o G/p " (8.3 X 10^ Ibo/hr-ft^)/ (0.179 Ibm/ft^) = 464,000 ft/hr = 128.8 
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The Reynolds Number is 
Re » Ug D^p/p 

» (464,000 ft/hr) (0.826 in.) (^^^j^) (0. 179 Ibm/f t^) /(0 .0889 Ibm/hr-ft) 
« 64,300. 

Using the circular tube friction factor correlation, 
f « 0.184/Re°-^ » 0.184/(64,300)°-^ 

- 0.0201. 

From the defining equation for the friction factor, 
APp = f(AL/D^) p Ug/2g^ 

= (0.0201)(21 ft)(r^—^)il20 in. /0. 826 in.)(0.179 Ibm/ft^) 

(464,000 ft/hr)^/(3600 sec/hr)^ (2) (32.2 Ibm-ft/sec^ - Ibf) 
= 282.8 Ibf/ft^ 
= 1.96 psi. 

1.96 

The average pressure in the core is 701 = 700 psia. Therefore, 

the assumed value of 700 psia was acceptable for t ^e density calculation. 
Example 6 

Each coolant channel in the reaCi.Ou core h^' the same pressure drop. 
If all coolant channels are not identic^^l t ie flew will redistribute such 
that this equal pressure drop is attained. Determine the relation 
between average velocity in the coolant channel (or equivalently the 

cf6' 
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mass velocity) and the diameter of 

distribution. 

Solution 

The frictional pressure drop 

AP « f 



the channel which governs the flow 



is given by 



Substituting for the friction factor. 



AP = 



0.184 AL 
ReO-2 2g^ 



0.184 



0,2 D 2g 
c °c 



D C 



if the same density and viscosity are assumed for each channel, 
channels 1 and 2, each having equal pressure drop. 



For 



-.1.8i_l,2 _-1,8j_1.2 
"bi /°c1 = "b2 



or 



"b2/"b1 = (°c2/°cl^ 



2/3 
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4.6 The Heat Transfer Coefficient 

No direct analytical methods for predicting the heat transfer 
coefficient for turbulent flow in circular tubes exist. However, by the 
use of an analogy between heat and momentum transport, the dependence of 
h on the flow field characteristics, i.e., on the Re can be predicted. 
Then by applying the results of an heat transfer analysis for flow over a 
flat plate an estimate of the Prandtl number dependence can be obtained. 
The Pr may be interpreted as a dimensionless modulus relating the temper- 
ature field in the fluid to the flow field. For Pr of unity, similar 
velocity and temperature profiles exist in the fluid. Finally, the 
predicted expression for h is compared to empirical correlations found 
to fit the available experimental data. 

The heat transfer across a fluid layer in laminar flow may be calcu- 
lated from Fourier^ s law 



where y is the direction normal to the fluid layer. Similarly, the shear 
stress in the fluid can be related to the velocity gradient by 



which defines the viscosity coefficient p. These relations can be 
rearranged to yield 



q = - k A 



(51) 



dy • 



T 




(52) 




(53) 



pCpA 



and 



T 

— = V 
P 



dU 
dy 



(54) 
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In these relations, a = k/pC^ is the molecular diffusivity of heat and v 
is the kinematic viscosity or the molecular diffusivity of momentum. We 
now postulate that the eddy transport of heat and momentum present in 
turbulent flow can be expressed in the same form. To this end, we define 
the eddy diffusivity of heat, e^, and the eddy diffusivity of momentum, e^^, 
as the parameters, that, when multiplied by the appropriate gradient, yield 
the corresponding transport rate for turbulent flow. In general, and 
throughout the flow field. For the combined molecular transport 
plus turbulent eddy transport, 

— ^ = O ^ (55) 



pCpA ^ '^H'' dy 

and 

We now assume that heat and momentum are transported by analogous 

processes and at the same rate. This requires that a ^ V and £„ = e^. 

H M 

Since the Pr = y/a, the first condition is equivalent to Pr « 1. Simi- 
larly a turbulent Prandtl number may be defined as Pr = ^ Therefore, 
the second condition is equivalent to Pr^ « 1. The basic assumption we 
have made also implies that both q and t vary in the same way across the 
flow field; i.e. , 

^ = Constant = — — — ^^^^ 



C At — — -CAT 
p p w w 



if the constant is expressed in terms of quantities at the surface or 
wall. 
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Dividing Equation 55 by Equation 56 and using the assumption dis- 
cussed above , 



CAT 
p w w 



dU « - dT, 



(58) 



Integrating this expression between wall conditions and fluid bulk 
conditions > 

iw_ r _ r (59) 



p w w U T 
*^ w w 



Using » 0, this reduces to 



= T - T 
CAT w B' 

p w w 



(60) 



To introduce h into this expression, Newton's law of cooling (Equation 27) 
is used to give 



A (T - T_) • 
W W B 



(61) 



Substituting Equation 61 into Equation 60, 



C T 

h - . (62) 

"b 



The shear stress at the wall can be eliiainated from this expression 
in favor of the friction factor. To do this, t^ is first related to the 
pressure loss over the length AL of the channel and then Equation 44 
relating the pressure loss and the friction factor is used. Consider 
the fluid in a. section of channel of diameter D and length AL. For fully 
developed flow, the sum of the forces acting in the flow direction must 
equal zero. The forces that must be considered are the pressure force 
on the upstream circular face, ^D^ P, the pressure force on the downstream 
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IT 2 ^ ■ 

face, (P-AP) , and the shear force along the periphery, it DAL t^. Sunmlng 

these forces with proper regard as to their direction gives 

y D^P = TT DALt + y (P-AP) (63) 
4 w 4 



or 



T =-5^AP 
w 4AL 



Using Equation 44 for AP, this becomes 



D . AL ^ B - ^ B f.,. 
w 4AL D 2 8 



Stibstituting Equation 64 into Equation 62, 

h = Cp p U3 I . (65) 

Introducing the Stanton Nrmber, a convenient dimensionless grouping 
common in heat transfer work. Equation 65 becomes 



P ^ 

Using Equation 49 for f, 

St = 0.023 Re''^'^. (67) 



This result is often found expressed in terms of the Nusselt nxjmber defined 
as 

Nu - P . (68) 

One common interpretation of the Nu is that it is the ratio of the actual 
convection heat transfer from a surface to the heat transfer assuming that 
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only molecular conduction were present. Noting that 

c. __h V! Nu_ (69) 

we can rewrite Equation 67 as 

Nu = 0.023 Pr Re^'^. (70) 

The Pr dependence found in Equation 70 is not correct. One reason 

why this incorrect dependence appears is that Pr = 1 was assumed early 

in the development; therefore, one would not expect the Pr dependence to 

be properly represented in the result. The analysis of laminar flow heat 

transfer from a flat plate, one of the few cases which can be treated 

analytically, yields a Pr"""^^ dependence is an expression for the Nusselt 

1/3 

number. In practice it turns out that this Pr dependence is reasonably 
accurate for the turbulent pipe flow case of interest here as well. In 
fact some of the empirical correlations found in the literature exhibit 
this Pr-"-^^ dependence in the Nu(Pr, Re) relation. However, the relation 
that is the best known, most widely used, and is recommended here is the 
Dittus-Boelter Equation [12] . 

Nu = 0.023 Re°-^ Pr". (71) 

For the fluid in the tube being heated the recommended value of n is 0.4 
and for the fluid in the tube being cooled the recommended value of n is 
0.3. The former case is the one found in reactor work. In Equation 71 
several temperature dependent fluid properties appear (k, p , y , C^) . 
These properties are to be evaluated at the average bulk temperature of 
the fluid for the segment of tubing of interest. The heat transfer 
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coefficient determined from Equation 71 is the average value around the 
circumference of the coolant channel. 
Example 7 

Determine the average heat transfer coefficient for the HTGR described 
in Exas:ples 1 and 5. 
Solution 

The average heat transfer coefficient corresponds closely to that 
calculated at the average coolant conditions: i.e., 1000 F and 700 psia. 
Additional fluid properties required for this calculation are 



C = 1.248 Btu/lbm-F , 



k =0.167 Btu/hr-ft-F. 



The Prandtl number is 



Pr = C m/V. 



= (1.248 Btu/lbm-F) (0.0889 Ibm/hr-f t) /0.167 Btu/hr-ft-F) 



= 0.664. 



From the circular tube Nusselt number correlation. 



Nu = 0.023 Re 



0.8 



0.4 



= 0.023 (64,300) 



0.8 



(0.664) 



0.4 



= 137.2. 



Calculation of h from the definition of Nu gives 



h = (Nu) k/Dj, 

= (137.2) (0.167 Btu/hr-ft-F)/ (0.826 in.)C 



12 in. 



= 332.8 Btu/hr-ft-F. 



ft. 
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4.7 Incremental Energy Balances 

As heat is transferred from the fuel to the coolant, the temperature 
of the coolant rises. Thus, the coolant temperature varies from a minimum 
at the inlet of the coolant channel to a maximum at the exit of the fore. 
The bulk temperature at any axial position can be determined from an energy 
balance where the total energy added by heat transfer to the coolant is 

equated to the energy rise of the coolant, i.e., 

Z 

il - ii, - r q' (Z) dZ (72) 
in / c 

■* H/2 

where q' is the heat transfer rate per unit length of fuel rod. Noting 
c 

that for a monatomic ideal gas Ai = C AT^, this can be rewritten as 

Z ^ 

h -hn-ld- /„,, "c "'^ 
p ^ H/2 

Equation 73 can be used to determine the bulk temperature increase for 

any segment, ^L, of the coolant channel by 

Z -AL/2 

1 . (74) 

T - T = / q'(Z) dZ 

B2 Bl m C 

c 

where Z is the center of the increment under consideration. Defining 
c 

an average surface heat flux for the increment, i^' = q^,/A^, the incremental 
bulk temperature rise becomes 

TT D AL 



T — T = r— a 

B2 Bl m C 

P 



c -.. (75) 



The local surface heat flux, q;'(Z), is related to the local volumetric 
thermal source strength by observing that all of the energy generated in 
the fuel is transferred into the coolant at the coolant-moderator surface, 

uD^ dZ q'^ (Z) = 2|d^ dZ q'"(Z). (76) 

5e 
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The factor of 2 on the right-hand side accounts for the fact that the 
energy generated in two fuel rods is transferred to each coolant channel. 
This simplifies to 

q^'(Z) - 2D ^'"(Z). (77) 



Similarly, 



c 



c 
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5.0 The HTGR Thermal-Hydraulics Code 

The HTGR Thermal-Hydraulics Code calculates the thermal-hydraulic 
performance parameters discussed in Section 4 for the HTGR model described 
in Section 3. The basic operational steps of the code are listed below. 

1. Real integer conversions 

2. Definition of statement functions 

3. Accept and print input 

4. Calculate geometrical parameters of reactor 

5. Calculate shape factors for finite-difference 
eqxiations used in moderator temperature 
distribution calculation 

6. Determine calculation increment 

7. Calculate inlet pressure loss. 

8. Initialize to 1st calculation increment 

9. Calculate bulk temperature rise and average 
temperature for 1st increment 

10. Calculate average coolant properties for 1st 
increment 

11. Calculate pressure loss, exit pressure, and 
average pressure for 1st increment 

12. Calculate heat transfer coefficient for 1st 
increment 

13. Calculate moderator temperature distribution for 
1st increment 

14. Calculate temperature distribution in fuel for 
1st increment 

15. Print 1st increment output 

16. Repeat steps 8-15 for remaining increments 

17. Calculate exit pressure loss and core exit pressure 

18. Print coolant exit conditions. 
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The details of these specific steps requiring further explanation 
are found in the discussion below. A flow chart showing the calculations 
and logic is given in Section 5.4. 

The statement functions are used for convenience in calculating 
property data. These statement functions are in the form of some property, 
either of the coolant, moderator, or fuel being the function of the 
corresponding temperature. All of the function forms are polynomials 
obtained by fitting polynomials of various degree to tabulated data. The 
order of these polynomials varies, having been selected to insure that it 
is accurate to written 2% over the range of interest of the parameter. 

A detailed discussion of the input requirements is deferred until 
Section 5.2. For the moment it suffices to note that this input is 
supplied on three cards; one each for the geometry, the inlet flow 
conditions, and the reactor power level. 

The temperature distribution in the moderator is calculated from 
the nodal equations developed in Example 4. As preliminary calculations, 
the coordinate locations of the nine nodal points and the conduction 
shape factors are determined from the coolant hole diameter, fuel diameter, 
and spacing supplied as input to the code. 

In the code, the model coolant channel, fuel rod, and moderator are 
sliced into short segments which are stacked axially to form the proper 
length core. The actual length of these segments is selected as a compro- 
mise between the very short segments for which the assumption of constant 
properties over the section axially is accurate and the very long segment 
which minimizes the computer time. An increment size equal to 1/50 of 
the active core length is built into the code. To study the effect of 



ERLC 



59 



56 



increment size, or to save computer time, this can readily be changed by 
altering two lines of the code. 

After the increment size has been computed, the coolant conditions at 
the beginning of the first increment are established. The bulk temperature 
and mass velocity at the entrance of the first increment are set equal to 
the core inlet conditions provided in the input. The pressure at the 
entrance of the first calculation increment is obtained by subtracting the 
inlet pressure loss, computed from Equation 44 and Equation 50, from the 
pressure at the inlet of the core. In addition, counters required in the 
code logic are set equal to 1, denoting the first increment. 

The first major calculation in the code is the determination of the 
bulk temperature rise for the first calculation increment. This temperature 
rise is calculated from Equation 75 with qjl determined from Equation 78. 
For small calculation increments, the average volumetric thermal source 
strength can be accurately approximated by evaluating Equation 1 at H/2 
-AL/2; i.e., at the center of the increment. 

The average bulk temperature, T_, for the increment is taken as the 

D 

arithmetic average of the inlet and exit values. All coolant properties, 
.except the density, used in the pressure drop and heat transfer coefficient 
calculations are assumed to be functions of temperature only and are 
evaluated at the average bulk temperature. The density is calculated from 
the ideal gas law using the average bulk temperature and the inlet pressure 
of the increment. 

Both the friction factor and heat transfer coefficient calculations 
require the Reynolds number of the flow which is calculated from Equation 
48. Then f and Nu are calculated from Equations 49 and 71, respectively. 




57 



Knowing f, the pressure drop for the increment is calculated from 
Equation 44. The heat transfer coefficient is calculated from the Nu 
definition. Equation 68. 

The temperature distribution in the moderator is calculated next. 
To do this, the nine finite-difference equations developed in Example 4 
are solved to yield the nine nodal temperatures. An elimination method 
is used in this calculation. The moderator temperature distribution 
calculation is complicated somewhat by the fact that both the moderator 
and fuel thermal conductivities which appear in the finite-difference 
relations are temperature dependent. Therefore, an iteration on the 
thermal conductivities is performed. The procedure is to assume average 
moderator and fuel temperatures, evaluate the thermal conductivities at 
these assumed temperatures from the statement functions, calculate the 
moderator and fuel (procedure described in next paragraph) temperature 
distributions and average temperatures, evaluate the thermal conductivities 
at these new average temperatures, and compare the new and original 
thermal conductivities. If the new thermal conductivities are within 
2% of those used in the previous iteration, the calculated moderator 
temperatures are accepted. If either thermal conductivity differs by 
more than 2%, the calculation process is repeated until finally the 2% 
criterion is met. The initial assumed average moderator and fuel tempera- 
ture are estimated from the results of the previous calculation increment. 
For the first increment, an initial guess is built into the code. 

After the coolant channel and moderator calculations for the incre- 
ment have been performed, the radial temperature distribution in the fuel 
is calculated for the increment. First, the temperature at the outer 
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surface of the fuel is found. This is taken as the area-weighted 

average of the temperatures calculated at the moderator-fuel interface 

in the moderator finite-difference calculation. The next step is to 

calculate the temperature at the center of the fuel. In this calculation 

axial symmetry is assumed and Equation 14 is used. In Equation 14 the 

thermal conductivity of the fuel appears. This property is temperature 

dependent and should be evaluated at the mean temperature of the fuel. 

However, the mean fuel temperature is dependent o\x the temperature rise 

across the fuel. Thus, an iterative-type of solution is called for. The 

procedure used is to evaluate at T^, calculate T^, calculate 

f , evaluate k at T , compare k to k . If 

fuel mean f-H fuel mean f-H ^ ^ 

k^ II is within 2% of J. accept value of T^. If k^_jj is not within 

07 of k let k = k and repeat the sequence. In this way, a 
°^ .'^f-i' f-i f-II 

fuel thermal conductivity accurate to within 2% is used. The fuel 
temperatures at the quarter, half, and three-quarter radii are then 
calculated from Equation 12. 

After the calculations for the first increment are completed the 
results are printed out. Discussion of the output is deferred until 
Section 5.3. It should be noted here, however, that in the output 
tabulation the temperatures are reported as having occurred at the axial 
center of the increment. In other words, for the first increment the 
results are reported at Z = H/2 - AL/2. After the results for the first 
increment have been printed out, the code increments to the next segment, 
and continues this process until all 50 segments have been spanned. 

After the completion of the calculations and printout for the 50 
segments of the core, the core exit conditions are calculated and printed. 
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A calculation of the exit pressure loss using Equations 44 and 50 is 

required. Then P and T are computed and printed, ending the compu- 
ex ex 

tatlons. 
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5 • 1 Nomenclature 



Analysis Symbol 
A 
A 

conv 



Code Symbol 



c^, c^, c^, 

COMPAl 
C0MPA2 



D 



e 
f 

G 
h 

H 



AFL 



C0MPA1 
C0MPA2 
CP 

DFU 
DCL 

ERR 
F 

GIN 
HTC 

H 



Description 

Area 

Convective heat transfer area 

Flow area of coolant channel 

Heat transfer area normal to n- 
directlon 

Heat transfer area at radius r 

Heat transfer area at surface 
of coolant channel 

Surface area of wall 

Integration constants 

Last iteration value of l^f(^^^) 

Last iteration value of k (T ) 

IE xn 

Constant pressure specific heat 
Diameter 

Outside diameter of nuclear fuel 
Diameter of coolant channel 
Energy release per fission 
Error ratio 
Friction factor 



Ibm-ft 



Dimensional constant - 32.174 2 



sec -Ibf 



Mass velocity 

Convective heat transfer 
coefficient 

Height of active core 



Si 



Analysis Symbol 
H 

e 

1 

INC 
k 

h 

k_ 



m 



'f-II 



8 

AL 

m 

n 

Nu 

Option 

Pw 
P 

P 

ex 
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Code Synbol 
HE 
ENTH 

INC 

KBC 

m0 

KFU 
KFUl 

KFU 2 



DELL 
MD0T 



NU 
0P 

PW 

PAV 

PI 

P2 

PEX 



Description 
Extrapolated height of core 
Enthalpy 

Enthalpy at inlet of core 

Increment counter in code 

Thermal conductivity 

Thermal conductivity of the 
coolant 

Thermal conductivity of moderator 

Thermal conductivity of fuel 

Thermal conductivity of fuel - at 
pre-iteration average fuel 
temperature 

Thermal conductivity of fuel - at 
updated average fuel temperature 

Flow resistance coefficient 

Extrapolation length 

Length increment 

Mass flow rate 

Exponent in Equation 71 

Number of fissionable fuel 
nuclei per unit volume 

Nusselt nuxnber 

Option = 1 for cosine. Option = 
2 for coupled program 

Wetted perimeter 

Pressure 

Pressure at beginning of increment 
Pressure at end of increment 
Core exit pressure 
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Analysis Symbol 
Pr 
AP 

q 

% 

q 

^gen 



Code Symbol 
PIN 
PR 



DELPE 
DELPF 



QTPAV 
QTFO 



Description 
Core inlet pressure 
Prandtl nuniber 
Turbulent Prandtl nunber 
Pressure drop 

Exit or entrance pressure loss 

Frictional pressure loss 

Heat transfer rate 

Heat transfer rate at 

Heat transfer rate in n-direction 

Heat transfer rate at D 

c 

Heat transfer rate at wall 

Heat transfer from node i to node j 

Heat flux at D 

c 

Energy generated inside nodal 
volume 

Volumetric thermal source strength 

Volumetric thermal source strength 
at center of core 



*old 
r 

r^ 



QTP(?LD 



RATI 

RAT2 

Re 
S 



RATI 

RAT2 

RE 
S 



Value of q"' used in previous 
increment 

Fuel rod radial coordinate 
Radias of fuel 

Core radial coordinate; also gas 
cons tant 

Relative error in fuel thermal 
conductivity between iterations 

Relative error in moderator thermal 
conductivity between iterations 

Reynolds number 

Center- to-center spacing between 
fuel rods and coolant holes 
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Analysis Synibol Code Symbol Description 

S. SI Coefficient of T. in nodal equations 

1 



Conduction shape factor between 
nodes j and k 



St Stanton number 

T T Temperature 

TO Temperature at center of fuel rod 

T^^^ T14 Temperature in fuel at r = 1/4 r^ 

T^^,j T12 Temperature in fuel at r = 1/2 r^ 

T^^^ T34 Temperature in fuel at r « 3/4 r^ 

T_ TBAV Coolant bulk temperature 

D 

Tg^ TBI Bulk temperature at beginning of 

increment 

Tg2 TB2 Bulk temperature at end of 

increment 

Tj TF Temperature in fuel at 

T^^ TF0AV Average temperature in fuel 

between rj/2 and r^ 

T TEX Core exit temperature 

ex 

T. TIN Core inlet temperature 

in *^ 

T TSAV Surface temperature inside convective 

layer 

T Free stream fluid temperature 

00 

T Wall temperature 
w 

AT Temperature difference 

U Velocity 

U UB Average velocity in flow channel 
B 

V Volume 

X Cartesian coordinate 

Coordinate in n-direction 
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Analysis Syinbol Code Symbol Description 

X Cartesian coordinate of nodal point 

AX Node spacing in x-direction 

y Coordinate normal to wall 

Y Cartesian coordinate of nodal point 
AY Node spacing in y-direction 

Z Axial coordinate 

Z ZC Axial location of center of incre- 

c 

ment 

a Thermal diffusivity 

e Roughness size 

Eddy diffusivity of heat 

Gj^j Eddy diffusivity of momentum 

0 Angular coordinate 

X^^ Transport mean free path 

y VIS Absolute viscosity 

V Kinematic viscosity 

p RH0 Density 

Average fission microscopic 
cross section 

T Shear stress; also time 

T Shear stress at wall 
w 

^ Neutron flux 

superscript- Average for calculation increment 
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5.2 Code Input 

The input to the HTGR Thermal-Hydraulics Code is intentionally 
extremely simple. Only three cards are required; the first lists the 
geometrical quantities, the second lists the coolant inlet conditions, 
and the third lists the parameters that describe the volumetric thermal 
source strength distribution in the reactor. 

The required input data cards are listed below. The units of each 
parameter and their input format are included. Sample input cards are 
shown in Figure 5. 

Card 1 - Geometry 

(in.) F 10.4 

(in.) F 10.4 

S (in.) F 10.4 

H (in.) F 10.4 

Card 2 - Coolant inlet conditions 

G (Ibm/hr-ft^) E 10.3 

T^^ (F) F 10.1 

P^^ (psia) F 10.1 

Card 3 - Power distribution 

q'" (Btu/hr-ft^) E 10.3 
o 

H (in.) F 10.1 

e 

Option (pure number) I 10 

On Card 3 the axial volumetric thermal source strength distribution 
through the core is specified. For independent operation of the code 
(Option = 1) a cosine distribution is built into the code and q'^' and 
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Figure 5. HTGR Thermal-Hydraulics Code Input. 
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He are the parameters which Indicate Its level and period. When the 
HTGR Thermal-Hydraullco Code Is coupled with other codes (Option « 2) 
the coslna distribution Is over-ridden by the actual power distribution 
which Is supplied by the main program. 
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5.3 Code Output 

Three sets of information are printed out. Each is. printed under 
completely explicit headings such that there should be no interpretation 
difficulties. 

The first set of output information consists of a listing of the 
input that was supplied to the code. 

The second and main set of output is a tabulation of the temperatures, 
h and AP for each increment in the core. The temperature printed out are 
f^, T^^^, T^^2» ^3/4» » ^B' '^1 ^^"^^ '^9 graphite moderator. 

From this data, observations of the maximum fuel and moderator temper- 
atures and their locations can be made. 

The third set of data reports the core exit pressure and temperature. 
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5.4 Cbdia Flow Chart 




Start 







kf(T), k^(T) 
kg(T), u(T) 







X 



Input: 
Dj, D^, S, H 



Functional relations between 
variables supplied to code as 
polynomials. 



Geometry, inlet coolant 
conditions > and power distribu- 
tion supplied to code. 



q'", H^, Option 



Output: 
Df, D^, S, H 

q^", H^, Option 



Printout of geometry. Inlet 
coolant conditions , and pover 
distribution. 



AL = H/50 



I 



1 = 1, 9 



Selection of calculation 
increment. 



Coordinates of nodal points 
for moderator finite-difference 
network are calculated. 
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jk = 12,14,15,23,25,26, 
36,37,45,56,58,67, 
68,69,79,89 



Output: 



Calculation of conduction shape 
factors for moderator finite- 
difference network. 



Printout of nodal point coor- 
dinates and shape factors. 



p = 


P. /RT. 
In In 




"b = 


= G/p 




Re = 




f = 


0.184/Re°-^ 




= 10f(AL/D^) Ug /2g^ 




= O.SpUg 





p = p 
1 In 

T T 
Bl = ^In 



^^E - ^^F 



Z = H/2 - AL/2 
c 

INC = i 



Calculation of pennament 
pressure loss at inlet. 



Initializing pressure, tempera- 
ture, position, and increment 
counter for .1st increment. 



T = + 300 cos ttZ /H^ 

m in c e 

= + 500 cos /H^ 

fo in c o. 



X 



Initial guesses for averL.;<« 
fuel snA moderator temperatures 
for first increment supplied 
to code. 

7.^ 
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, • t • 



q'" cosirZ /H 



^B2 - Tg^+q'"irD^/2fiC 



Calculation of average bulk 
temperature of coolant for 
increment. 







" G/p 




Re » U_D p/ 

o C 




f = 0.184/Re • 


APp - f(AL/D^)pUg/2g^ 


P« «■ P - AP 
2 1 F 



Pr » l,248y(Tg)/k2(Tg) 
Ku = 0.023Pr°-^Re°-^ 
h - Nu lcB(Tg)/D^ 



1 VSi<^^,i.H<^fo>'W'^'^>j 



Calculation of pressure drop 
and exit pressure for increment. 



Calculation of ^average heat 
transfer coefficient for 
increment. 



Obtaining best guess initial 
values for average fuel and 
moderates ten^eratures for 
finite-difference calculation. 



Cjilculation of coefficients cl- 
T,. in nodal equations. 
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q'^.Dj.D^.h.Tg) 



» (T^+Tg) /2 



1 - <^f) 



Solution of moderator teinperature 
distribution. 



Calculation of average temperatures 
at surfaces of fuel and coolant. 



Calculation of fuel centerllne 
teiq>erature. The thermal con- 
ductivity of the fuel Is eval- 
tiated at a mean fuel teo^erature 
by an Iterative technique. 



T J - (Tj+T^) /2 



- k,(T^) 




COMPAl - kj(Tj^) 

C0MPA2 » k (T J 
m m 



Saving values of thermal con- 
ductivities for later comparison. 
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T^^-(3T^+5Tp/8 

^m" l'*"^2''"^3''"^4''"^5''"^6'*"^ y"*"^ 8^"^ 9 ^ 

RATI - |(COMPAl-kj(Tj^))/COMPAl| 

RAT2 » |(C0MPA2-k (T ))/C0MPA2l 
mm ' 




fl 



^3/4'^o-5«l' "4/256 k,^ 



Output: 

1* 2' 3' 4' 5' 
Tg.T^.Tg.Tg 




Check to deterndne if correct 
moderator and fuel thennal 
conductivities were used in 
moderator calculation. If not, 
an iteration on the average 
thermal conductivities is 
performed. 



Evaluation of fuel ten^erature 
at the quarter, half, and three 
qtiarter diameter points. 



Printout for increment. 



Check to see if entire core 
has been calculated. 



78 



74 



INC - INC+1 



Z ■ Z -AL 

c c 

T = T « 
Bl B2 




Initializing for next Increment. 



Ug-G/p 

Re-UgD^e/p(Tg2) 
f=0.184/Re°*^ 

F 



AP„-10f(AL/D^)pUg/2g^ 



APg-pU3/2g^ 



- 


^2- 


•^^E-^^F 


T 

ex 


T 

^B2 



Output: 



P , T 
ex ex 



Calculation of exit pressure 
losses. 



Calculation of coolant core 
exit conditions. 



Printout of coolant core exit 
conditions • 
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5.5 Code Examples 

The following two examples demonstrate the use of the HTGR Thermal- 
Hydraulics Code for calculating the thermal-hydraulic behavior of a 
typical HTGR. 

Example 8 . 

A HTGR core is composed of 3000 right hexagonal prism fuel elements, 
each 16 in. across the flats and 32 in. high. The elements are stacked 
eight high in the core. The uranium- thorium fuel in carbide form fits 
tightly into 0.8 in. diameter fuel holes. The reactor is cooled by the 
downward flow of helium which enters at 500 F, 600 psia and flows through 
1.0 in. diameter coolant holes. The spacing between the fuel and coolant 
holes is 1.25 in. The average coolant velocity in the inlet of the core 
is 280 ft/sec. The average fuel rod produces 4.2 kw/ft. Assimie an 
extrapolation length equal to 5% of the core height. Determine the 
magnitudes and locations of the maximum fuel temperature and the maximum 
moderator temperature. 



The mass velocity, extrapolated core height, and peak volumetric 
thermal source strength must first be computed from the data given. 



Solution. 



(600 ^)(144^) 



G - p, U. 



U, 




in ft 



(280 — )(3600 



sec 



) 



in B in R T 



B in 



-) 1.500 + 460) R 



sec 



hr 



in 



» 2,35 X lO"^ Ibm/hr-ft 



H -H+2L -K+ 2(0. 05H) - I.IOH 



e e 
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= (1.10)8(32 in.) = 281.6 in. = 23.47 ft 

q't. = 2P /[H sin(TrH/2 H )] 
o ave e f ^ e' ^ 

= 2(4.2 kw/ft)(8)(32 in.) ( ) (3413 Btu/kw-hr) /{ (23 .47 ft)(1.0 in.)^ 

12 in. 

(ft/12 in.)^ sin(Tr(8)(32 in.)/2(28l.6 in.)]} 
= 3.79 X 10 Btu/hr-ft"'. 



From the output printout of the code (next 7 pages) we find, 



T = 889.7 P at Z = -3.63ft at node 1 

mod max 



. = 962.9 F at Z = -2.77ft. 
fuel max 
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HTGR THERMAL-HYDRAULICS CODE 
INPUT DATA 



FUtL 0.0 COOL O.D. 

(INCH) (INCH) 
0.8000 1.0000 



SPACING 
(INCH) 
1.2500 



MASb VtLOCITY 
(L8/HR-FT**2) 
0.2350E C6 

VCLU>tETRIC THERMAL 
SOURCE STRENGTH 
(BTU/HR-FT**3) 
0.3790E C7 



CCRE INLET TEMPERATURE 
(F) 
500.00 

EXTRAPOLATED HEIGHT 

( INCH) 
281.60 



ACTIVE CORE LENGTH 
(INCH) 
256.00 

CORE INLET PRESSURE 
(PSIA) 
600.00 



OPTION 1 



t N D tJF INPUT DATA 



DFU=0.G667 



COORDINATES OF THE NODES 

DC0=0.0833 S=0.1042 



NJOL NUMBER 

1 

2 
3 
4 
5 
6 
7 
3 
9 



X COORDINATE 

0.18750E-01 
0.23216E-01 
0.35417E-01 
O.OOOOOE 00 
0.11608E-01 
0.11608E-01 
0.20833E-01 
O.OOOOOE 00 
O.OOOOOE 00 



Y COORDINATE 

O.OOOOOE 00 
0.16667E-01 
0.28868E-01 
O.OOOOOE 00 
0.16667E-01 
0.28868E-01 
0.54127E-01 
0.28868E-01 
0.48544E-01 
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P A T ».4 


SHAPh FACTlJk 


o H D 


J. 5625C0CB 


jO 




0.31i7972E 


00 




0.4d221<.CE 


00 




0.2142fc28fc: 


00 


S12h 


0.866025^£ 


00 




0.116101&E 


01 






00 




0.^757009E 


00 




0.8660254E 


00 




0.9514C17E 


00 




O.A75/009E 


00 


^ . .. V 


VJ.2949650E 


00 


Sod 


0. i37i0y8L 


01 


S Jo 


0.4310531E 


00 


So7 


0.427692CE^ 


00 


S6y 


0.3509365E 


00 


SJ7 


0. 182ol62t 


00 


S79 


0. 4106029b 


JO 


AT IHl 


iNLLT dF THE CURb 





CHA;\i><EL AVE. VELOCITY KEYNOLDS # FRICTION FACTOR 
(FT/SEC) 

279.97 280285.9 O.Olt.0 
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NU Z-LUC TO Tl/4 Tl/ii T3/4 TF 

TN(I) TN(2) TN(3» TNt^l TN(5) TN(6» 
TN(7I TN(a) TN(9) TSURF TBULK HTCOF OELP 

1 10.453 551.7 551.2 548.2 543.8 537.7 

538.7 537.7 536.5 536.9 535.8 533.1 

523.7 532.6 525.4 525.6 500.6 690.9 21.86 

2 10.027 569.8 569.1 565.1 559.3 551.1 
552.5 551.2 549.6 550.1 548.7 545.1 

533.3 544.5 535.0 535.1 502.0 691.0 21.90 

3 9.600 588.4 587.5 582.5 575.2 565.. 0 

566.8 365.1 563.1 563.7 562.0 557.4 
543.3 536.7 544.8 545.0 503.8 691.1 21.96 

4 9.173 607.1 606.1 600.1 591.4 579.2 
581.3 579.3 576.8 577.6 575.6 570.1 

f,.,j.4 569.2 '554.8 555.1 506.0 691.2 22.02 

b 8.747 626.0 624.6 617.9 607.8 593.6 

396.0 593.7 590.8 591.7 589.3 582.9 

563.8 561.9 565.1 565.5 508.5 691.4 22.09 

O t.320 643.0 643.6 635.7 624.2 608.1 

610.9 608.2 605.0 606.0 603.2 595.9 
376.3 594.7 575.6 376.0 511.5 691.6 22,17 

/ 7.893 663.9 662.4 653.6 640.7 622.7 

623.9 622.9 619.2 620.3 ^17. 3 609.0 
387.0 607.7 586.2 586.6 514.8 691.9 22.27 

Q 7..467 682.8 681.2 671.5 657.3 637.4 

O40.9 637.6 633.5 634.8 631.3 622.2 
597. o 620.7 596.9 597.4 518.4 692.1 22.37 

9 7.040 701.6 699.8 689.2 673.8 652.1 

633.9 652.3 647.8 649.2 645.5 635.4 
6U8.7 633.8 607.7 608.2 522.4 692.4 22.48 

10 6.013 720.2 718.3 706.9 690.1 666.7 
670.9 667.0 662.1 663.6 «>59.5 646.6 

619.6 646.9 618.5 619.1 326.7 692.7 22.60 

11 6.187 738.6 736.5 724.3 706,4 681.3 
683.8 681.5 676.3 677.9 ^'^3.5 661.8 

63U.6 659.9 629.3 629.9 331.3 69:».0 22.72 
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1^ ■ b./60 7b6.o 75^.^ 741,4 72^,5 095.7 

7oj,b . o9:>»9 690.4 692. J 6ti7.4 674.8 

041.4 0/2.8 6-^>o«l 640.8 i;>36.2 6v3.3 22«io 

IJ :j.iJ3 774.3 771.9 75fc.2 738.0 709.9 

7lH.y 710.2 704.2 706.0 Vul.l 687.7 
6b^.3 6b5.6 650. d <>51.5 b41.4 693.7 23.00 

14 4.907 791,3 789.0 774. < 753.4 723.8 
7Z9.2 724.1 717.9 719. d 714.5 700.4 

o63.U 69d.2 661.5 662.2 546.9 694.1 23.15 

15 4.4d0 803.1 805.6 790.5 768.4 737.5 
743. 1 737.8 731.2 ?33.2 727.7 712.9 

073.5 710.6 672.0 67i.7 552.6 o94.5 23.31 

iu H.053 824.2 82i,6 805. S 782.9 750. d 

750.0 751.1 7-V4.3 746.3 740. 6 725.1 
ub3.9 722.6 682.0 683ci 55b. 5 694.9 2j.4/ 

17 3.o27 839.7 836.9 820-7 796.9 763.7 

7o9.d 7o4.0 756.9 759.0 753.0 736.9 
09'!. 1 734.4 692.4 693.2 564.7 695.3 23.64 

lb 3.200 854.5 851.6 ii34.'; dlG.4 776.2 

7d2.5 776.5 769.1 771.3 765.1 748.4 
7w4.0 745.8 702.3 703.2 571.0 695.7 23.81 

U 2.772 863.5 865.6 (>48.4 323.3 783.1 

794.6 788.5 780.9 783.1 776.7 759.5 

713.7 750. d 711.9 712.8 577.5 696.1 23.99 

20 2.347 881.7 878.7 b6l.2 ' 835.5 799.5 
dU6.2 799.9 792.1 794.4 787.8 770.2 

723.1 767.4 721.2 722.2 584.1 696.6 24.17 

21 1.920 894.1 891.0 673.1 847.0 81C.t 
ai7.2 810.8 802.8 805.1 798.4 780.3 

7:)2. 1 777.4 730.2 731.2 590.9 697.0 24.36 

22 1.493 9C5.5 902.4 884.3 857.7 820.6 
827.5 821.0 812.9 815.2 tOB.4 789.9 

740.8 787.0 73<J.8 739.8 597.8 697.5 24.55 

23 1.067 ~l6.u 912.9 894.5 867.7 830.1 
8J 7.2 830»':' 822.3 824.7 817.7 799.0 

749.1 796.0 7^7.1 748.1 604./ 697.9 24.74 
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24 0.640 925.5 922.4 903.9 876.8 839.0 

846.1 839.4 831.1 833.4 826.4 807.5 
756.9 804.4 754.9 755.9 611.7 698.4 24.94 

25 0.213 934.0 930.8 912.3 885.1 847.1 
854.3 847.5 839.1 841.5 834.4 815.3 

164.4 812.3 762.4 763.4 618.8 698.9 2*^.13 

26 -0.213 941.4 938.3 919.7 892.5 854.5 
861.6 854.9 846.5 848.8 841.7 822.5 

771.3 819.5 769.3 770.3 625.8 699.3 25.33 

27 -0.640 947.8 944.6 926.1 899.0 861.0 

863.2 861.4 853.0 855.4 848.3 829.0 
777.8 826.0 775.8 776.8 632.9 699.8 25.52 

2b -1.067 953.0 949.9 931.5 904.5 866.8 

873.9 867.2 858.8 861.1 854.1 834.9 
783. ii 831.8 781.8 782.8 639.9 700.2 25.72 

29 -1.493 957.1 954.0 935.8 909.0 871.7 
878.8 872.1 863.7 866.1 859.0 840.0 

7b9.2 837.0 787.2 788.2 646.8 700.7 25. 9i 

30 -1.920 960.1 957.0 939.0 912.6 875.7 
882.8 876.1 867.9 870.2 863.2 844.4 

794.1 841.4 792.2 793.1 653.7 701.1 26.10 

31 -2.347 961.9 958.9 941.2 915.? 878.9 
885.8 879.3 871.2 873.5 366.6 848.0 

798.5 845.1 796.5 797.5 660.5 701. t 26. 'j9 

32 -2.773 962.5 959.6 942.^ 916.8 881.2 

888.0 881.6 873.6 875.9 869.1 £50.9 
302.3 848.0 800.3 801.3 667.1 702.0 26.48 

33 -3.200 962.0 959.2 942.2 917.4 882.6 

889.3 883.0 875.2 877.4 870.8 853.0 
805.5 850.1 803.6 804.5 673.6 702.4 26.66 

34 -3.627 960.4 957*6 941.1 917.0 883.2 
889.6 363.6 876«0 878.? 871.7 854.3 

80a.l 851.6 806.3 807.2 680.0 702.8 26.84 

35 -4.053 957.6 954.9 938.9 915.5 882.8 

889.1 883.2 875.9 877.9 871.7 854.9 
810.1 852.2 808.3 809.2 686.1 703.2 27.01 
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3o -4,t3u 953.7 951.0 935.7 Sil3.1 881.6 

887.6 882.0 674.9 876.6 670.9 854.7 
811.5 852.1 a09.8 810.6 692. U 703.6 27.13 

37 -4.907 948.6 946.1 931.3 909.7 879.5 

885.3 879.8 873.1 874.9 869.2 853.7 

812.3 851.2 810.6 811.5 697.7 703.9 27.34 

33 -5.333 942.5 940.1 926.0 905.4 876.6 

882.1 876.9 870.4 872.2 866.7 851.9 

012.5 849.6 010.9 811.7 703.2 704.3 27.50 

39 -5.7O0 935.3 933.0 919.7 9U0.i 872.8 

878.0 873.1 866.9 868.6 863.4 849.4 

812.0 847.2 810.5 811.3 708.4 704.6 27. > 

^0 -6.187 927.1 924.9 912.3 893.9 866.1 

873.1 868.4 862.7 864.3 859.4 846.2 

811. U 844. 1 809.6 810.3 713.3 704.9 27.79 

^1 -6.613 917. b 915.8 904.0 886.8 862.7 

867.3 663.0 857.6 859.1 854.5 842.2 
809.3 840.2 80a. 0 80b. 7 717.9 705.2 27.93 

42 -7.040 907.6 905.8 894.8 878.9 856.5 

600.7 656.7 851.7 853.1 848.9 637.5 
6u7.i 835.7 805.9 806.5 722.2 705.4 28.05 

^3 -7.467 896. D 894.6 884.8 870.1 849.5 

853.4 849.7 845. 2 846.4 842-6 832.1 

804.2 830.4 803.1 603.7 726.2 705.7 28.17 

44 -7.893 884.5 883. C 873.9 860.5 841.8 
845.4 842.0 837.9 839.0 835.5 82^*-,. I 

800. 3 824.5 799.8 800.3 729.9 705.9 2ii»2ci 

45 -8.320 671.7 870.3 662.2 0.2 833.4 
830.6 833.6 829.9 831.0 827. 819.4 

796.8 818.0 795.9 796.3 733.2 .06.1 28.33 

46 -8.747 858.1 856.9 649.7 339.2 824.4 

827.2 824.6 821.3 622.2 819. J; 812.0 
792.2 810.8 791-4 791.8 736.1 706.3 28.47 

47 -9.173 843.8 842.6 830.6 827.5 814.7 
817.1 814.9 812.1 812.9 t;10.5 804^1 

787.1 8C3.1 786.4 766.7 738.7 706.4 28.55 
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48 -9.600 828.8 8rj,0 822.8 815.1 804.5 

806.5 804.6 802.3 802.9 801.0 795.6 
781.4 794.8 780.9 781»l 740.8 706.6 ZU.bl 

49 -10.027 813.3 812.5 808.4 802.3 793.7 
795.3 793.8 791.9 792.5 790.9 786.6 

775.3 785.9 774.8 775.0 742.6 706.7 28.67 

50 -10.453 797.1 796-6 793.4 7fe8.8 782. V 

783.6 782.5 781.1 781.5 780.3 777.1 
766,6 776.6 768.3 768,4 744.0 706.8 28.72 



EXIT PRESSURE EXIT TEMPERATURE 

(PSIAJ <FJ 

584.2 744.6 
POWER PLR UNIT LfNGTK (KW/FTJ s 2.69 
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Example 9 , 

The effect of undersized coolsnt holes on HTGR thermal^hydraullc 
performance is to be studied. Consider the HTGR described in Example 8 
but with a coolant hole diameter of 0.8 in. (instead of 1.0 in.). 
Compare the following thennal-hydraulic parameters for the "undersized 
coolant hole" and "normal" flov? channels: 

a. Location and magnitude of maximum fuel temperature. 

b. Location and magnitude of maximum moderator temperature. 

c. Exit coolant ttiiapetature . 

Solution . 

Each flow channel in the core will experience the same pressure 
drop. The "undersized coolant hole" channel has a smaller diameter and 
hence will pass less flow for the same pressure loss (as shown in 
Example 6) . To assure approximately the same pressure loss as in the 
"normal" channels, the mass velocity of the "undersized coolant hole" 
channel miat be adjusted. From the results of Example 6, 

2/3 

%ch^S%ch/V 

= 2.35 X 10^ (0.8/1.0)^/^ 
= 2.025 X 10^ Ibm/hr-ft^ 

The inpi.t to the HTGR Thermal-Hydraulics Code is the same as that 
Cor Example 8 except for the new value of G and the undersized D^. 
From the code (printout next 7 pages) the following results are ob- 
tained: 

a. The maximum fuel temperature is now 1189 F and occurs at 
Z « -3.63 ft compared with 963 F at Z = -2.77ft for the 
"normal" channel. 
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HTGR THERMAL-HYDRAULICS CODE 
INPUT DATA 



FUEL O.D COOL O.D. 

(INCH) (INCH) 
0.8000 0.8000 



SPACING 
(INCH) 
1.2500 



MASS VELOCITY 
(LB/HR-FT**2) 
O.2025E C6 

VOLUMETRIC THERMAL 
SOURCE STRENGTH 
(BTU/HR-FT**3) 
0.3790E 07 



CORE INLET TEMPERATUKE 
(F) 
500.00 

EXTRAPOLATED HEIGHT 

( INCH) 
281.60 



ACTIVE CORE LENGTH 
(INCH) 
256.00 

CORE INLET PRESSURE 
(PSIA) 
600.00 



OPTION 



END OF INPUT DATA 



Of U=0.0667 



COORDINATES OF THE NODES 

DCO=0.0667 5=0.1042 



NODE NUMBER 

I 

2 
3 

5 
6 
7 
8 
9 



X COORDINATE 

0.18750E-01 
0.23216E-01 
0.35417E-01 
O.OOCOCE 00 
0.11608E-01 
0.1160aE-01 
0.16667E-01 
O.OOOOOE 00 
O.OOOOOE 00 



Y COORDINATE 

O.OOOOOE 00 
0.16667E-01 
0.20868E-01 
O.OOOOOE 00 
0.16667E-01 
0.28868E-01 
0.61343E-01 
0.28868E-01 
0.56878E-01 
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SHAPE FACTORS. 



PATH 


SHAPE FACTOR 


S45 


a.5625CO0E 


00 


S14 


0.3117972E 


00 


S15 


0.4822120E 


00 


S12M 


0. 2142628E 


00 


5J.2F 


0.8660234E 


00 


S25 


0. lL6i010E 


01 


S26 


0.4999998E 


00 


S23M 


0.4757009E 


00 


S23F 


0«8660254E 


00 


S56 


0.9514017E 


00 


St>8 


0.4757009E 


00 


sa9 


0.2072C95E 


00 


S66 


0.1732048E 


01 


536 


0.6331320E 


00 


S67 


0.3545S11E 


00 


S69 


0.2210881E 


00 


S37 


0.7788467E- 


-01 


S79 


0.8202949E 


00 


AT THE 


INLET UF THE CORE 





CHANNEL AVE. VELOCITY PEYNLLDS U FRICTION FACTOR 
(FT/StC) 

241.25 193218.4 O.OlOl 
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NO Z-LOC TO Tl/4 Ti/2 T3/4 TF 

TNfl) TN(2I TN(3J TNC4) TN(5) TN(6l 
TN(7» TN(8J TN(9) TSURF TBULK HTCOF DELP 



1 10.453 564.7 564.2 561.2 556.8 550.7 

551.7 550.7 549.5 549.8 548.8 546.0 
534.8 545.6 534.6 534.7 501.1 641.3 21.86 

2 10.027 588.3 587.6 583.6 577.8 569.6 
570.9 569.6 568.1 568.5 567.1 563.3 

548.4 562.8 548.1 5^*6.3 503.6 641.5 21.93 

3 9.600 612.6 611.7 606.7 599.4 589.2 
590.9 589.2 587.3 587.8 586.0 581.3 

562.5 580.6 562.2 562.4 506.9 641.7 22.02 

4 9.173 637.4 636.4 630.4 621.6 609.3 
611.4 609.4 607.0 607.6 605.5 599.9 

577.2 599.0 576.8 577.0 510.8 642.0 22.13 

5 8.747 662.6 661,4 654.4 644.2 629.9 
632.4 630.0 627.3 628.0 625.5 618.9 

692.3 617.9 591.8 592.1 515.5 642.3 22.26 

6 8.320 688.2 686.8 678.9 667.3 651.0 

653.8 651.1 647.9 648.7 645.9 638.3 

607.8 637.2 607.3 607.5 520.8 642.6 22.40 

7 7.893 714.0 712.5 703.6 690.6 672.4 
675.6 672.5 669.0 669.9 666.7 658.1 

623.7 656,8 623.1 623.4 526.8 643.0 22.57 

B 7.467 740.0 738.4 728.6 714.2 694.1 

697.6 094.3 690.3 691.3 687.8 678,2 

639.9 676.7 639.2 639.5 533.4 643.4 22.75 

9 7.040 766.1 764.3 753.6 737.9 716.0 

719.9 716.2 711.8 712.9 709.0 698.5 
656.3 696.9 655.6 655.9 540.6 643.9 22.94 

10 6.613 792.2 790.3 778.7 761.8 738.1 
742.3 738.2 733.5 734.6 730.4 719.0 

672.9 717.2 672.2 672.5 548.4 644.4 23.15 

11 6.187 818.2 816.1 803.7 785.6 760.1 

764.7 760.4 755.2 756.4 751.9 739.5 
689.7 737.6 688.8 689.3 556.8 644.9 23.38 
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12 5>.760 3A<»eO 841.8 309.3 782.2 
787.1 782.4 776. S 77d.2 773.3 760.0 

700. i. 758.0 705.6 706.1 bOt>.7 645.5 23.62 

13 5.333 869.5 867.2 853.2 832.8 604.2 

809.4 804.4 798.6 799-9 794.7 780.5 
723.4 778.3 722.4 722.9 575.1 646.0 23.88 

14 4.907 894.7 892.2 b77.5 856.1 826.0 

831.5 826.3 820.1 821.5 816-0 800.9 

740.3 798.6 739.3 739.8 5b5.0 646.7 24.14 

15 4.480 919.4 916.7 901.4 879.0 847.6 

853.3 847.9 841.3 842.8 637.0 821.1 

757.0 818.6 756.0 756.5 5<*5.A 647.3 24.43 

16 4.053 943.4 940.7 924. b 901.5 868.8 
b74.b 8o9.1 8t>2.2 863.8 o57.6 840.9 

^75.7 838.4 772.6 773.1 6u6.1 648.0 24.72 

17 3.627 966.9 964.1 947.6 923.4 889.6 
895.9 tt89.9 882.8 884.3 677.9 860.5 

790.1 €57.8 789.0 789.5 617.3 648.6 25.02 

10 3.200 989.6 986.7 969.6 944.8 909.9 

916.4 910.2 902.6 904.4 897.7 879.6 

806.4 876.8 805.2 805.6 628.7 649.3 25.33 

19 2.773 1011.4 1008.4 990.9 965.4 929.6 
936.4 929.9 922.3 923.9 917.0 898.2 

822.3 895.3 821.0 821.6 640.5 650.0 25.65 

20 2.347 1032.3 1029.2 1011. A 985.3 948.7 

955.6 949.0 941.1 942.8 93i).7 916.2 
837.9 913.2 836.6 837.2 652.6 650.8 25.98 

21 1.920 1052.2 1049.1 10J0.9 1004.3 967.0 
974.1 967.4 959.2 96J.9 953.6 933.7 

tt53.1 930.6 851.8 852.4 6o4.6 651.5 26.32 

22 1.493 1C71.0 1067.8 10'*9.4 1022.4 934.6 
991.8 984.9 976.6 978.3 970-8 950.4 

b67.8 947.2 966.5 867.2 677.3 t>52.2 26.66 

23 1.067 1088.6 1085.4 lOftfc.7 1039.4 1001.2 
lOOtf.6 1J01.6 993.1 994.8 987.2 966.3 

S82.1 963.1 880.8 881.4 639.9 653.0 27.01 
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24 0.640 1105.0 1101.8 1083.0 1055.5 1016.9 

1024.4 1017.3 1008.7 1010.4 1002.7 981.5 
895.9 978.2 894.5 895.2 702.6 653.7 27.36 

25 0.213 1120.1 1116.9 1098.0 1070.3 1031.6 
1039.2 1032.0 1023.4 1025.1 1017.2 995.8 

909.1 992.4 907.7 908.4 715.4 654.4 27.71 

26 -C.213 1133.8 1130.6 1111/7 1084.0 1045.3 
1052.9 1045.6 1036.9 1038.6 1030.7 1009.1 

921.7 10C5.7 920.2 921.0 728.2 655.2 28.06 

27 -0.640 1146.1 1142.9 1124.0 1096.4 1057.8 

1065.5 1058.2 1049.4 1051.1 1043.2 1021.5 
933.6 1018.1 932.2 932.9 741.0 655.9 28.42 

28 -1.067 1156.9 1153.7 1135.0 1107.5 1069.1 

1076.8 1069.5 1060.8 1062.4 1054.5 1032.8 

944.8 1029.4 943.4 944.1 753.7 656.6 28.77 

29 -1.493 1166.3 1163.1 1144.5 1117.3 1079.3 

1086.9 1079.6 1071.0 1072.6 1064.7 1043.0 
955.4 1039.6 953.9 954.7 766.3 657.3 29.12 

30 -1.920 1174.0 1170.9 1152.6 1125.7 1088.1 
1095.7 1088.5 1079.9 1081-5 1073.7 1052.1 

965.1 1048.8 963.7 964.4. 778.7 658.0 29.47 

31 -2.347 1180.2 1177.2 1159.1 11>32.7 1095.7 

1103.2 1096.1 1087.6 1089.2 1081.4 1060.2 
974.1 1056.8 972.7 973.4 791.0 658.7 29.81 

32 -2.773 1184,8 1181.8 116^.1 1138.2 1102.0 

1109.3 1102.3 1094.0 1095.5 108719 1067.0 
982.3 1063.7 980.9 981.6 803.1 659.4 30.15 

33 -3.200 1187.8 1184.8 1167.6 1142.3 1106.9 
1114.1 1107.3 1099.1 1100.6 1093.1 1072.6 

989.6 1069.4 988.3 989.0 814.8 660.0 30.48 

34 -3.627 1189.2 1186.3 1169.5 1144.9 1110.5 
1117.5 1110.9 1102.9 1104.3 109''. 0 1077.1 

996.1 1073.9 994.8 995.5 826.3 ^60.6 30.81 

3b -4.053 1188.9 1186.1 1169.9 1146.1 1112.8 

1119.5 1113.1 1105.3 1106.7 1099.7 1080.3 
1001.7 1077.2 1000.4 1001,1 837.5 661.2 31.12 
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3o -4.4U0 1167.0 1184. A 1168.7 1145.7 1113.6 

1120.2 1113.9 1106.5 1107.8 1101.0 1082. J 
lu06.4 1079.3 1005.2 1005.8 848.2 661.8 31.43 

37 -A. 907 1183.6 1181.0 1166. C 1143.9 1113.1 
1119,4 U13r.4 1106.3 1107.6 lltJl.O 1083.0 

1010.2 1080.2 1009.0 1009.6 858.6 662.3 31.73 

38 -5.333 1178.5 1176.1 1161.7 1140.7 1111.3 

1117.3 1111.6 1104.8 1106.0 1099.8 1082.6 
lOlj.l 1C79.9 1012.0 1012.5 668.5 662.9 32.01 

39 -5.7)0 1171.9 1169.6 1156.0 1136.1 1100.2 
1113.8 1108.4 1102.0 1103.1 1097.2 1080.9 

1015.0 1078.4 1013.9 1014.5 677.9 663.4 32.28 

40 -6.187 1163.8 1161.6 1148.8 1130.0 1103.7 
11C9.1 1104.0 1097.9 1099.0 1093-4 107tt.l 

1016.0 1075.7 1015.0 1015.5 866.8 663.8 32.54 

41 -6.613 1154.3 1152.2 1140.2 1122.6 1098.0 
11J3.0 1098.2 1092.5 1C93.6 1088.3 1074.0 

1016.1 1071.8 1015.1 1015.6 895.2 664.3 32.78 

42 -7.040 1143.2 1141.3 1130.2 1113.9 1091. J 
1095.7 1091.3 1086.0 1086.9 1082.1 1068.8 

1015.2 1066.8 1014.3- 1014.7 903.0 664.7 33.01 

43 -7.467 1130.9 1129.1 1118.9 1103.9 1082.9 
1087.1 1083.1 1078.2 1079.1 1074.7 1062.5 

1013.3 1060.6 1012.5 1012.9 910.2 665.0 33.22 

44 -7.893 1117.2 1115.6 1106.3 1092.6 1073.5 

1077. 4 1073.7 1069.3 1070.1 1066.1 1055.1 
1010.6 1053.4 1009.8 1010.2 916.8 665.4 33.42 

4i> -8.320 1102.2 1100.8 1092.4 1080.2 1063.1 

1066.5 10b3.3 1059.3 1060.1 1056.5 1046.6 
1006.9 1045.1 1006.2 1006.5 922.3 665.7 33.59 

46 -8.747 1086.1 1084.8 1077.5 1066.7 1051.6 

1054.6 1051.7 1048.3 1048.9 1C45.8 1037.2 
1002.3 1035.8 1001.7 1002.0 928.1 666.0 33.75 

47 -9.173 1066.6 1067.7 1061.4 1052.1 1039.1 

1041.7 1039.2 1036.2 103u.8 1034.1 1026.7 
996.7 1025.5 996.2 996.5 932.7 666.2 33.90 
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4a -9.600 1050.5 1049.6 1044,3 1036.5 1025.6 

1027.7 1025.7 1023.2 1023.7 1021.4 1015.3 
990.J 1014.3 989.9 990.1 936.7 666.4 34.02 

49 -10^027 1031.2 1030.5 1026.2 1019.9 1011.2 
1012.9 1011.3 1009.3 1009.7 1007.9 1003.0 

s;83.i 1002.2 982.7 982.9 939.9 666.6 34.12 

50 -10.453 1011.0 1010.5 1007.2 1002.5 996.0 
997.2 996.0 994.5 994.8 993.5 989.8 

974.9 989.2 974.7 974.8 942.5 666.7 34.21 



EXIT PRESSURE EXIT TEMPERATURE 

(PSIAJ (F) 

583.5 943.6 



POWER PER UNIT LENGTH (KW/FT) : 2.69 
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b. The maximum moderator temperature is now 1120 F and occurs at 

node 1 at Z = -4,48ft compared with 890 F for node 1 at Z = -3.63 
ft for the "normal" channel. 

e. The helium exit temperature is now 944F compared to 745 F 
for the "normal" channel. 
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5.6 Listing of Code 



c 

C VPI HTGR THERMAL-HYDRAULICS CODE. 

I 

c 

c 

KtAL KFU,KM0,KdCtN0fMD0T,KFUl,K|-U2 
DIMENSION T(9)fA(81)iB(9} 

C 

|C TtMPERATUKE DEPENDENT PftOPEHTlES. 

C 

KFU(T)=0.1fi483E*02-0.72l35E-02*I+0.46776E-05*T«*2-C.10256E-08* 
< T»*3 

KMO(T)=0.53309E+02-0.31556E-QI*T+0.6629E-05*T**2 
K6C(T)=0.8ll425E-01+0.100714E-03*T-0.1A2851E-07*T**2 
VISIT J=u,40903E-OUO,6274E-04*T-0.96155E-OB*T**2 

C READING INPUT DATA. 

ftEA0(5,5)DfUiDC0,S,H 
5 FORMAT(3Fl0.4,Fi0.n 

R£AO(5,10)GIN,TIN,PIN 
10 FQRMATI£10.3,2F10.1) 

R£A0(5,15)UTP0,HE.0P 
15 FaRHAT(El0.3fF10.l.I10) 

I PRINTOUT OF INPUT DATA. 

HRITEl6,20) 

20 FGRMAT(lHlf2X///////i5Xf'**«»***************************/15X, 
<«HTGR THERMAL-HYDRAULICS C0DE'/15X, • 



<«//i5X, 'INPUT DATA'//) 
WRITE(6,25) DFU,DCO,SiH 
25 FORMAT {ISXt'FUcL fJ.DS6X, 'CluL CO. ' , 7X, • SPACINd' i 7X, • ACTI VE CORE 
< LENGTH* ,6X,/15X,M INCH) ' ,9X , • ( INCH)S9X, 
<• tlNCHIM4X,« nNCH}M4X/15X,F6.4,9X,F6.4,9XtF6.4t 
<14XfF6.2tl6X/) 
WRlT£(6,30)GIN,Tn-;,PIN 
30 F0KMAT(15X,'«ASS VEIOCI TYS6X, 'CORE INLET TEMPERATURE • ,3X, 
<«CJRE INLET PR£SSURfc«/15X,'(LLN/fiK-FT=?*2)M3X,' (F) 't20X,' IPSIA)'/ 
<15XtElU4,14X,F7.2,l7X,F7.2/) 
WRITE|6,35)QTP0,HE 
35 F0RMAT(15Xt»V0LUHETRIC THERMAL ' ,4X, 'EXTRAPOLATED HEIGHT • /15Xt ' SOUR 
<Ce STR£NGTH»/15X,MBTU/HR-FT'p*3)MlX,MlNCHI»/l5X,tU.4,14X, 

<F6.2//) 
kyRlTE(6,40)OP 

40 F0RHAT(l9Xf •GPTI0NS3X,I2///15X, 't N D OF INPUT DATA',/ 
<15X|- »,////) 

C 

C CJNVERTING INCHES TO FEET* 

DFU=DFU/U.O 

S=S/12.0 

DCU=0C0/L2.0 

H=H/i2.0 

hE=HE/12.0 

C 

C A SYSTEM UF NODES IS btT UP FOR THE FINITF-DIFFEKENCt ANALYSIS 

C OF HEAT TRANSFER IN THE GRAPHITE. THE ORIGIN OF THIS SYSTEM, 

C THE CENTER OF THE FUEL, AM; ThE CENTER OF THE COOLANT HOLE FORM 

C A RIGHT-ANGLED TRIANGLE. THE 'XSANO 'Y' COORDINATES OF THE 

C NODES ARE CALCULATED IN TERMS OF S,THE SPACING, OFU,THE DiA- 

C METER OF THE FUEL, AND l,)CO,THt DIAMETER OF THE COOLANT HOLE 
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C AS FOLLOWS: 

C 

PI=ARC0S(-1.0) 

Xl=S*C0S(Pl/3.0)-0,5*DFU 

Y1=0.0 

X2=S*C0S(Pl/3.0)-0«5*DFU*C0S(PI/6.0) 

Y2=0.5WU«SINIPI/6.0) 

X3=S*COS(Pl/3.0)-0,5*OFU*COS( PI/3.0) 

Y3=0.5WU*SIN(PI/3.0) 

X4=0.0 

Y4=0.0 

X5=X2/2.0 

Y5«Y2 

X6=X5 

Y6=Y3 

X7*0.5*0CC*$IN(PI/6.0) 

Y7=S*COS(PI/6.0)-0.5*OCO*COS| PI/6.0) 

X8:=0.0 

Y8=Y3 

X9=0.0 

Y9=S*CGS(PI/6.0)-0.5*DC0 

C 

C THE FOLLOWING SEGMENT CALCULATES THE CONDUCTION SHAPE FACTORS 

C ASSOCIATED WITH THE HEAT TRANSFER BETWEEN THE NODES. EACH SHAPE 

C FACTOR CORRESPONDS TO HALF THE SUM OF THE COTANGENTS OF THE 

C ANGLES ENCLOSING THE PATHS. 

C 

S45«0.5*(X5/Y5*(X1-X5)/Y5) 

S14=0.5*(C0TANIATAN(X5/Y5)+ATANnxl-X5)/Y5)J) 

S15=0.5*(X5/Y5+IX2-Xl)/Y2) 

S12M=0.5*(XI-X5)/Y5 

S12F=0.5*C0TAN(PI/6.0) 
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i)23=U.5*!CGTAN(ATAN({X2-Xli/Y^i)+ATAM(Xl-X51/Y51l+COTAK(ATWa(X2 

< -X6}/(Y6-Y2)H-ATANItXb-X5l/{Y6-Y5)))) 
S26=0.5*(IX6-X3)/IY6-Y5)+(X3-X2)/(Y3-Y2)) 
S23M=C.5«I(X2-X6)/(Y6-Y2)) 
S2iF=0.5*COTANlPI/6.0) 
S56=G.5*(IX2-X6)/|Y6-Y2) + Xi)/(Yd-Y5)) 
S58=0.5*(tX6-X5)/iYb-Y5) +X5/(Y8-Y3)) 
S39=0.5*X6/(Y9-Y6) 

St)3=0•5*((Y9-Y8)/X6+CaTA.\(ATA^:(X5/(Y3-Y5ll+ATAN{IX^)-Xp)/ 

< (Y6-Y5)))J 

S36=0.5»(C0TAiM(ATAN((X7-X6)/(Y7-Y6l)tATAN((X3-X7)/(Y7-Y3))J+ 

< C0TAN(ATANl(X2-X6)/{Y6-Y2)l+ATAMIX3-X2)/(Y3-Y2))l) 
$67=u*5*{CDTAN(PI-ATAN{X6/(y9-Y6))-ATAMX7/iY7-Y9)))+ (X3-X7)/ 

< (Y7-Y3I) 

So9=u.5«CCTAN(PI-PI/3.0-ATAi\'{(Y7-V9)/X7)-ATAN(lX3-X7)/IY7-Y31)- 

< ATAN((X7-X6)/IY7-Y6)II 
S37=0.5*(X7-X6)/IY7-Y6) 

S79=0.5«CCTANIPI-ATAN((Y9-Yfc)/X6l-ATAN((Y7-Y6)/(X7-X6))) 

C 

C PRINTOUT OF THE CClORDINATtS AND THE SHAPE FACTORS. 

C 

WRlTE(6i45) OFUiDCOtS 
'^j FORMAT(/30X,'CCUKOINATES OF THE ,U'DlSV/15X,' DFU=' tF6•'^^, 
<10Xf*DCC=',Fo.^,lOX,'S=SF6.4//) 

WRITE 16,30) Xl,Yl,X2,Y2,X3,Yi,X4,Y4tX5,Y5,Xo,Y6,X7,Y7,X8,Y6,X9,Y9 
5G FORMAT115X,'MtjOE NUMBER" ,3Xi ' X CUUROINATE* ,9Xt*Y COOROINATEV/^OX, 

< •lMOX,Ell.5tl3X,E11.5/ZOX,'2MJX,E11.5,10X,GiU5/20X,'3'tia, 

< Ell.5,10X,tii.5/2OX,»4MLX,EU.tjil0X,Ell.5/2OX,»5M0X,Ell.5,nx 

< ,Ell,5/20X,'oM0X,Ell.Si0XTElU5/20X,'7M0X,EU.5iL0X,Ell.W 

< 2oX,'8M0X,hU.5,l0Xtt;ll.^/2GX,»9M0X,ElU5,10X,Ell.5/) 

wRIT£(6,55) 
b'j FQRMATliHl,2X//////) 
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WRITE (6,60) S45,S14,S15,S12rt,S12F,S25fS26,S23M,S23F,S56,S58,S89, 

< S68,S36,S67,S69,S37,S79 

60 FORMAT (22X, •SHAPE FACTORS. V/Z/ISX, 'PATHMOX, 'SHAPE FACTOR*// 

< 15X,»S45M0X,E14,7/15X,»S14M0X,£14.7/15X,«S15M0X,EH.7/ 

< 15X,»S12M», 9X,£14.7/15X,»S12F»t 9X»E14.7/15X,»S25M0X,E14.7/ 

< 15X,«S26M0X,E14.7/15X,'S23MS 9X,E14.7/15X,«S23F' , 9X,E14.7 

< /15X,*S56M0X,£14.7/15X,»S58M0X,E14.7/15X,»S89M0X,£14.7/ 

< 15X,«S68MOX,E14.7/15X,«S3eM0X,El4.7/l5X,«S67M0X,£14.7/ 

< 15X,»S69M0X,E14.7/15X,»S37M0X,E14.7/15X,»S79M0X,E14.7///I 

C MASS FLOW RATE CALCULATION. 

MOOT=GIN*PI*OCO**2/4.0 
C SELECTION OF AXIAL CALCULATION INCREMENT. 

0ELL=Hi^'50.0 

C CALCULATION OF PRESSURE LOSS AT INLET OF CORE. 

RHOsPIN*144.0/(TIN«460.0)/3d6.0 
UB=6IN/RHQ/3600.0 
RE=UB*OCO*RH0/VIS(TIN»*3600.0 
Fs0.184/RE**0.2 

DELPF=F*DELL*RH0*UB**2/0C0/2. 0/32.2*10.0 
DELPE«0.5*RH0*UB**2/2.0/32.2 

C 

C PRINTOUT OF CORE INLET CONDITIONS. 

C 

WRITE(6,65) 

65 FORMAT! 15X, 'AT THE INLET OF THE CORE*//) 
WR1TE(6,70) UBiREiF 



1^3 



N FURMAT{iDX,»ChANN£L AVE. VELLXITY • ,2X, 'REYNaiDS ,3X, 'FRICTION ^ 
<ACTOR'/21X,'(FT/SECP/21X,F6.2,I2X,Fc1.1t8X,F7.'V//I 

KRITE(6,75) 
<S///) 

c 

L PRINTOUT CiF COLUMN htAi)I(vioS. 

C 

3j FORMAT(IHh2X///////l5X,'(\0'f7X,'Z-LGC»t5X,»TO't6X,'Tl/4»,5X, 
<'Tl/2»,5X,«T3/4' ,7X,«TF'/20X,»TrJ(l)»,4X,»TN(2)S4X,'TNl3P,4X, 
<»TN{4)»,4X,«TN(5JS4X,«TN(6)«/l5Xt»TN(7)*f4X,«TN(8)S4Xi»TNI9)», 
< 5XfnS0RF',4X,'TBULKS4X,»HTCOFS5X,»DELP»//) 
^RITE(6i8C) 

C 

C CALCULATION OF CONDITIONS AT hEGINNING CF FIRST INCREMENT, 

C 

Pl=PIN-iO£LPE+OELPF)/l44.0 

TS1=TIN 

INC=1 

C 

t CALCULATIuN OF AXIAL COOPUINATE UF CENTER OF FIRST INCREMENT. 

C 

ZC- H/2.0-0tLL/2.0 

C 

C INITIAL GUESS OF AVERAGE FUEL AND MCOERATOK TEMPERATURES FOR 

I THE FIRST INCREMENT. THESE ARE USED TO EVALUATE THE THERMAL- 

C CONDUCTIVITIES IN THE CALCULATION OF THE MODERATOR TEMPERATURE 

C DISTRIBUTION. 

C 

TMAV=TIN+30C.0*CLjSIPI*ZC/rit) 
TFOAV=TIN+50u.0*COSlPI*ZC/HE) 
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c 

C CALCULATION OF AVERAGG BULK TEMPERATURE FOR INCREMENT. 

C 

85 OTPAV=QTPC*CO$(PI*ZC/HEl 

TB2=TBl+QTPAV*Pl«DFU**2.0*DELL/4.0/M00T/U24800*2,0 

C 

C THE MULTIPLICATION BY 2 IN THE PREVIOUS CALCULATION IS REQUIRED 

C BECAUSE TWICE THE ENERGY GENERATED IN ONE FUEL ROD ENTERS EACH 

C COOLANT CHANNEL. 

C 

TBAV=(TBl+TB2)/2.0 

C 

C CALCULATION OF PRESSURE DROP AND EXIT PRESSURE FOR INCREMENT. 

C 

RHa=Pl*144.0/ITBAVH60.0)/386.0 
UB=GIK/RHO/3600.0 
RE=UB*OCO»RH0/VIS(TBAV)*3600.0 
F«0.184/RE««0.2 

OELPF=F*DELL*RHO*UB**2/0CO/2.0/32.2 
P2=P1-0ELPF/144.0 

C 

C CALCULATION OF HEAT TRANSFER COEFFICIENT FOR INCREMENT. 

C 

PR=1.248*VIS(TBAV)/KBC(TBAV) 

NU=0.023*RE**0.8*PR**0.4 

HTC=NU*KBClTBAV)/DCO 

C 

C CALCULATION OF TEMPERATURE DISTRIBUTION IN GRAPHITE MODERATOR. 

C 

C 

C 

C OBTAINING THE BEST GUESSES FOR AVERAGE FUEL AND MODERATOR 
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C TEHPtRAmES FOR Fl^aT::-DI^FhKENCt: CALCULATION JF MODERATOR 



C 



r 



TEMPERATURE DISTRIBUT luN. 



IF(INC.EQ.l) GO TG 95 

TMAV=TBAV*(TMAV-T6AV)*CTPAV/ijTPr.LC 
TFJAV^TBAV+(TF0AV-TtiA\/|-7^TPAV/QTPCLD 
95 RATIO=KF'J(TFaA\/l/ltMC(TKAV) 

C 

C CALCUUTINl, THE MAIN uIAbui^AL COEFFICIENTS FOR THE COEFFICIENT 

C MATKIX UF THE NCDAL EQUATK](\S. 

i 

S1=S1A+S15+S12M+RATI0*S12F 

S2=S26+S25+Si2M+S12F*RATia+S23M+RATlG*S23F 

S3=S36+S37+S23M+RATI0*S23F 

S5=S45+S15*S25+S56*S58 
S6=S56+S26+S36+S67+S69+S63 
S7=S37+S67+S79+PI*0CQ*hTC/{24.0*KMG(TMAVI) 
S8=Si)8+S68+S09 

S9=S89*S69+S79+{PI»0C0*HTC/(24.C*kMC(TMAV))) 



C SETTING UP THE CLtFFIClCi\T >.ATRIX FOR SIMQ SOLUTION OF NCDAL- 

C EUIATIONS. 

C 

00 100 I=lf9 

100 Bin=o.o 
00 1C5 i=ifei 

103 AII)=0.0 
A(l)=-Sl 

AliO)=lS12M+RATIb*S12r) 

. A(28)=SU 
A(37)=S15 



1^6 



ERIC 



A(2)=A(10) 
A(ll)=-S2 

A(20)=(S23M+RAT10*S23F) 

A(38)=S25 

A(47)=S26 

A(12)=A(20) 

A(48I=S36 

A(57}=S37 

A(21)=-S3 

A(4}=S14 

A(40)=S45 

A(31)=-S4 

A(5)=S15 

A(14)»S25 

A(50)-S56 

A(66)-S58 

A|4l)=-S5 

A(32)=S45 

A(42)=S56 

A(15)=S26 

A(24|5S36 

A(60)=S67 

Aa8)»$69 

A(69)=S68 

A(51)=-S6 

AI25)=S37 

A(52)=S67 

A(79)=S79 

A(61l=-S7 

A(71)=-S8 

A(44)*S58 

A(53l=S6d 



A(80)=S89 
A(81)=-S9 
A(72)=Sd9 
A(d4I=S69 
A(63I=S79 

b(U=-QTPAV*HI*Ll'U*'^2/(9u.C-KJ'.Li(rf<AV)) 
fi<3)=BIU 

8(7)=-PI*0Ca?HrC*TBAV/(24.0*^Kf''ClTMAVI) 
B(9)=B(7I 

KKS=0 

C 

C USING Sim TO OBTAIN SOLUTION OF NODAL EQUATIONS. 

C 

CALL SIMQ(A,BtNNNfKKS) 
DO 110 1=1,9 

110 T(n=B(n 

c 

C CALCULATION OF AVERAGE TEMPERATURE AT SURFACES OF FUEL AND 

C COOLANT FOR INCREMENT. 

C 

TF=(Ta)+2.0'^T|2)*T(3))/4.j 
TSAV=lTl7)+T(^))/2.0 

L CALCULATION OF TEMPEKATUKb AT CENTER Or FUEL FOR INCRtMLNT. 

C ITERATION TECHNIQUt b USlD TC OBTAIN FUEL THERMAL CONDUCTIVITY 

C WITHIN 2«. 

C 

KFJl=KFUiTF) 
115 TU=TF+gTPAV*LlFiJ*«2.0/16.J/KFUi 
TFAVS={TFnC)/2.0 
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KFU2=KFU(TfAVS) 
ERR-ABS({KFUl-KfU2)/KFU2) 
IF(ERR-0.02 )125tl25,12C 
.120 KFU1=KFU2 
Gu TO 115 

« 

C 

C SAVING VALUES OF THERMAL CONDUCTIVITIES FOR LATER COMPARISONS. 

C 

125 C0MPA1=KFU(TF0AV) 
C0MPA2=KMC(TMAV) 

C 

C CHECK TO DETERMINE IF CORRECT MODERATOR AND FUEL THERMAL- 

C CONDUCTIVITIES WERE USED IN MODERATOR FINITE-DIFFERENCE 

C CALCULATIONS. IF NOTtAN ITERATION ON THE AVERAGE THERMAL- 

C CONDUCTIVITIES IS PERFORMED. 

C 

TFQAV=3.0/8.0*TO+5.0/8.0nF 

TMAV=(T(l)fT(2l*T(3)+Tl4)+T(5)+T(6)+T(7)*T(8)fTI9))/9.0 
RATl=ABSIiC0MPAl-KFU(TFOAV))/C0MPAll 
RAT2=ABS( (C0MPA2-KH0(TMAV) )/C0MPA2l 
IF(RAT1«LE. 0.02. AND. RAT2.LE. 0.02) GO TO 130 

GO TO 95 

C 

C CALCULATION OF FUEL TEMPERATURE AT 1/4, 1/2, AND 3/4 FUEL RADII 

C FOR INCREMENT. 

C 

130 T14=T0-QTPAV*DFU»*2.2/16.0/KFU1/16.0 
T12=T0-QTPAV*DFU**2/16.0/KFUl/4.0 
T34=TC-QTPAV*DFU**2. 0/16. O/KFUl/16. 0*9,0 

C 
C 

C PRINTOUT OF RESULT FOR CALCULATION INCREMENT. 
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c 

WRITE 16,135) INC,ZC,TOJl4tT12,T34,TP,B,TSAV,TBAV,riTC,0ELPF 
U5 F0RMAT(15X,I2WX»F7.3,2X»F6,i,4X,F6.1,3X,F6.l,3X,P6.1,4X,F6.1/ 
<20X,F6.l,5(3X,F6.l)/15X,3IF6.l,3X),lX,F6.1,2(3X,F6.U,4X,F6,2/) 
iF(INC.EQ.ll.0R.INC.eQ.23.0R.INC.EQ.35,0R.INC,EC1.47) WRITE(6,55I 

C 

C CHECK TO SEE IF CALCULATIONS HAVE BEEN COMPLETED FOR ENTIRE 

C "CORE. 

C 

IFIlNC-50) 14Cfl45fl45 

C 

C SETTING CONDITIONS FOR BEGINNING OF NEXT CALCULATION INCREMENT. 

C 

140 INC=INC+l 
ZC=ZC-OELL 
TB1=TB2 
Pl=P2 

QTPOLD=CTPAV 

C 

C TRANSFERRING TO START CALCULATIONS FOR NEXT INCREMENT. 

C 

60 TO 85 

C 

C CALCULATION OF CORE EXIT CONnillONS OF COOLANT. 

L 

145 RH0=P2*l44.C/(TB2+460. 01/386.0 

UB=GIN/RHO/3600.0 
R£=UB*OCO«kHO/VIS(TB2)*3600.0 

F=0.iB4/RE**0.2 

DELPF=F*DELL*RHU*UB**2/DCO/2.0/32.2«lC.O 

DELPE=RH0*Ue**2/2.0/32.2 
PEX=P2-(DELP£+DELPF)/144.0 
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TEX=T62 

C 

C PRINTOUT OF CORE EXIT COOLANT CONDITIONS. 

C 

WRITEl6i75) 
MRITE(6iI50) 

150 F0RMAT(/iI5X,»EXIT PRESSUREMOX,»EXIT TEMPERATURE' ,/19X, 
<MPSIA)M8X,» (F)»/) 
MR1TE(6,I55)PEX,TEX 
135 F0RMATa8X|F6.Ifl8X|F6.I) 

C 

C CALCULATION OF AVERAGE POWER OF FUEL ROD IN KW/FT. 

C PRINTOUT OF AVERAGE POWER. 

C 

P0VL=0TP0*H£WU**2/2.0»SINIPI*H/2.0/HE)/H 
P0VL»P0VL/3412 
URITE(6«160)POVL 
160 F0RMAT(//15X, 'POWER PER UNIT LENGTH (KW/FT) :SF6.2/////) 
URITE(6J5) 
WRITE(6i55) 
STOP 
END 

C 

C 

C 

C SUBROUTINE SIHO 

C 

C PURPOSE 

C OBTAIN SOLUTION OF A SET OF SIMULTANEOUS LINEAR EQUATIONSi 

C AX=B 

C 

C USAGE 
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CALL SIMQ(A|biN,KS) 



DESCRIPTION OF PARAMETERS 

A - MATRIX OF COEFFlCIttiTS STORED COLUMNWISE. THESE ARE 
DESTROYED IN THE COMPUTATION. THE SUE OF MATRIX A IS 
N BY N. 

B - VECTOR OF ORIGINAL CONSTANTS (LENGTH N), THESE ARE 

REPLACED bY FINAL SOLUTION VALUESi VECTOR X, 
N - NUMBER OF EQUATIuNS AND VARIABLES. N MUST BE .GT. ONE. 
KS - OUTPUT DIGIT 

0 FOR A NORMAL SOLUTION 

1 FOR A SINGULAR SET OF EQUATIONS 

REMARKS 

MATRIX A MUST BE GENERAL. 

IF MATRIX IS SINGULAR , SOLUTION VALUES ARE i»IEAN INGLE SS. 
AN ALTERNATIVE SOLUTION MAY BE OBTAINED BY USING MATRIX 
INVERSION' (MINV) AND MATRIX PRODUCT (GMPRD). 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 

METHOD 

METHOD OF SOLUTION IS BY ELIMINATION USING LARGEST PIVOTAL 
DIVISOR. EACH STAGE OF ELIMIf>3ATICN CONSISTS OF INTERCHANGING 
ROWS WHEN NECESS/iKY TO AVCIO DIVISION BY ZERO OR SMALL 
ELEMENTS. 

THE FORWARD SOLUTION TO ObTAIN VARIABLE N IS DONE IN 
N STAGES. THE BACK SOLUTION FOR THE OTHER VARIABLES IS 
CALCULATED BY SUCCESSIVE SUBSTITUTIONS. FINAL SOLUTION 
VALUtS A^E DEVELOPED Il\ VICTOR b» WITH VARIABLE I IN BlDi 
VARIABLE 2 IN B(2)t t VARIABLE N IN BIN). 
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C IF NO PIVOT CAN 6E HOUND EXCetOING A TULfcRANCE OF 0.0» 

C THE MATRIX IS CONSIDEftED SINGULAR AND KS IS SET TO 1. THIS 

C TOLERANCE CAN BE MODIFIED BY REPLACING THE FIRST STATEMENT. 

C 

C 

c 

SUBROUTINE SIMQl A,B|N,KS) 
DIMENSION A(6I),6(9) 

C 

C FORWARD SOLUTION 

C 

T01.=0.0 

KS=0 

JJ=-N 

DO 65 J=1»N 

JY=J+1 

JJ=JJiNfl 

BIGA=0 

IT=JJ-J 

DO 30 U^fH 

C 

C SEARCH FOR MAXIMUM COEFFICIENT IN COLUMN 

C 

IJ=IT*I 

IF(ABS(BIGAJ-A6S(A(IJ))) 20i30i30 
20 BIGA=A(IJ) 

IMAX=I 
30 CONTINUE 

C 

C TEST FOR PIVOT LESS THAN TOLERANCE (SINGULAR MATRIX) 

C 

IF(ABS(BIGA)-TOL} 35>3S,40 



35 KS=1 
RETUkN 

INTERCHANGE ROWS IF NECESSARY 

40 I1=J+N*(J-2I 
I T= I MAX- J 
00 50 K=J,N 
Ii=Il+N 
I2=I1+IT 
SAVE=A(I1) 
A(il)=A(I2) 
A{12)=SAVE 

DIVIDE EQUATION BY LEADING COEFFICIENT 

50 A(IU=A(I1I/BIGA 
SAVc^BdMAX) 
B(IHAX)«B(J) 
B(J)==SAVE/BIGA 

ELIMINATE NEXT VARIABLE 

IF(J-N) 35,70,55 
55 IQS=N«U-1) 

DO 65 IX=JY,N 

IXJ=ICS+IX 

1T=J-IX 

DO 60 JX=JY,N 

IXJX=N*(JX-I)+IX 

JJX=IXJX+IT 
60 A(lXJX)=A(IXJX)-IA(IXJi*^A(JJX)) 
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o5 6iIX)=B(IX|-(b(J)«A(IXJ)J 

BACK SCLUTION 

70 NY-N-1 
IT=N*N 

00 80 J=lfNY 

IB=N-J 
IC=N 

DO 80 K=1,J 

e(IB)=fin6)-A(IA)«B(IC) 

80 IC=IC-1 
RETURN 
END 
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C 
C 
C 
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Problems 



By performing an energy balance on a differential element with sides 
dx, dy, and dz, derive the steady state heat conduction equation in 
Cartesian coordinates. 

Derive the steady state heat conduction equation in cylindrical 
coordinates . 

Develop an expression for the temperature distribution in a long, 
thin, hollow cylindrical fuel rod with inner radius r^, outer radius 
r^, and volumetric thermal source strength q'". Assume that the 
inner surface is insulated and at temperature T^. 

A long, thin, hollow cylindrical fuel rod with volumetric thermal 
source strength q'" is cooled on both the inner and outer surfaces 
(r - r^ and r - r^ surfaces) such that the surface temperatures are 
T. and T , respectively. Derive an expression for the location of 
the maximum temperature in the fuel. 

Calculate the ratio of peak power to average power for a fuel rod of 
constant cross section assuming that the volumetric thermal source 
strength varies as q'" » q'l' cos (TrZ/H^) and that the extrapolation 
length is 10% of the core height. Compare this ratio to that for the 
case where the extrapolation length is assumed to be zero. 

Develop an expression for the steady state heat transfer rate through 
a cylindrical shell of inside radius r^ and outside radius r with 
no internal energy generation. Your expression should contain only 
the outside surface temperature T , the inside surface temperature 
T^, the thermal conductivity k, the length of the shell AL, r^, and 

A HTGR fuel rod is 0.6 in. in diameter and 30 in. long. The surface 
temperature of the rod is at 1200 F and the volumetric thermal source 
strength of the rod is 8 x 10^ Btu/hr-ff^. The thermal conductivity 
is 12.0 Btu/hr-ft-F. Determine the temperature and heat transfer 
rate at r = 0.2 in. 

The fuel rod in Problem 7 fits inside a fuel hole that is 0.630 in. 
in diameter. The 0.030 in. diametrical gap between the fuel rod and 
moderator is filled with helium gas. Assuming that the heat transfer 
through the helium is by molecular conduction only (k^^ - O'l^O 
Btu/hr-ft-F) determine the temperature difference between the fuel 
and moderator surfaces. 

Develop the steady state finite-difference nodal equation for node 2 
of Figure 2. 

Develop the steady state finite-difference nodal equation for node 4 
of Figure 2. 
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11. Develop the steady state finite-difference nodal equation for node 10 
of Figure 2. 

12. Develop the steady state finite-difference nodal equation for node E 
of Figure 3. 

13. Develop the steady state finite-difference nodal equation for node H 
of Figure 3. 

14. Develop the steady state finite-difference nodal equation for node J 
of Figure 3. 

15. A 25 foot long HTGR core contains 0.520 in, diameter coolant holes. 
The average conditions in the core are 500 psia and 1000 F, The mass 
velocity of the helium is 2.0 x 10^ Ibm/hr-ft^. Determine the pres- 
sure drop through the core. 

16 • Determine the effect on pressure drop of a 20% decrease in coolant 
hole diameter for the core in Problem 15. Assume the coolant condi^ 
tions and mass velocity remain the same. 

17. Determine the heat transfer coefficient for the flow conditions of 
Problem 15. 

18. A HTGR core contains 0.75 in. diameter coolant holes through which 
helium flows. At a particular location in the core the pressure is 
500 psia, the temperature is 900 F and the coolant velocity is 216 
ft/sec. Determine the heat transfer coefficient. 

19. In a HTGR core the volumetric thermal source strength varies 

axially as q''' (Z) = q ' ' ' cos (irZ/H ). The energy generated in 

each two fuel^rods is removed by a coolant channel which has a mass 

flow rate of m and an inlet temperature of T^^. Designate the fuel 

diameter as D^ and the coolant hole diameter as D . Develop an 

expression for the axial variation of the temperature at the surface 

of the moderator. Your expression should contain only T^^, q^''', 

D., D , H, H , m, C , h, and Z. ^ 
f c e' p 

20. From the expression developed in Problem 19 determine the axial lo- 
cation where the moderator surface temperature (moderator-coolant 
interface) is a maximum. 

# 

21. A 4000 Mw (t) HTGR contains 5400 right hexagonal prism fuel elements 
each 14.2 in. across the flats and 31.2 in. long. The fuel elements 
are stacked 8 high in the core. The uranium- thoriiin fuel* in carbide 
form is contained in 0.62 in. diameter fuel holes. The core is 
cooled by the downward flow of helium through 0.825 in. diameter 
coolant holes. The helium enters^at 600 F, ^25 psia, and with an 
average mass velocity of 8.3 x 10 Ibm/hr-ft . Each fuel element 
assembly contains 132 fuel holes. The spacing between the fuel and 
coolant holes in the element is 1.112 in. The thermal source 
strength varies as q'" « q coa ^Z/H and H « I.IH. Determine 
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the maximum fuel temperature and the maximum moderator temperature 
in the core, 

22. Determine the % increase in power that would produce a 1800 F 
maximum fuel temperature for the core and flow conditions of Problem 
21. 

23. Determine the % reduction in core flow that would produce a 1800 F 
maximum fuel temperature for the core and flow conditions of Problem 
21. 

24. A 2000 Mw (t) HTGR contains 2400 right hexagonal prism fuel assemblies 
stacked 10 elements high to form the core. Each fuel element is 

12 in. across the flats and 30 in. high. In each element the uranium- 
thorium fuel in carbide form is contained in 150 fuel holes each 0.48 
in. in diameter. The core is cooled by the downward flow of helium 
through 0.56 in. coolant holes. The helium enters at 700 F, 600 psia, 
and with an average velocity of 300 ft/sec. The spacing between the 
fuel and coolant holes is 0.80 in. The axial power distribution is 
a trtjncated cosine ftmction with an extrapolation length of 16 in. 
Determine the location and magnitude of the maximum fuel, moderator, 
and coolant temperatures. 

25. The coolant holes in one of the fuel element assembly columns of the 
core in Problem 24 are undersized by 30%. Determine the % increase 
in maximum fuel, moderator, and coolant temperatures over those for 
the normal reactor fuel element assemblies. 

26. Determine the % overpower that leads to a maximxjm fuel temperature 
of 1800 F for the core and flow conditions of Problem 24* 
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